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ABSTRACT 


A  two-dimensional  finite  element  model  of  a  low  frequency, 
free- flooded  piezoelectric  thin  ring  sonar  transducer  was 
developed  for  use  with  the  ATILA  code. 

Effective  material  properties  for  the  ring  were  determined 
from  in- air  modal  analyses  using  ATILA.  These  were  adjusted  to 
obtain  the  closest  agreement  between  the  calculated  ring 
resonance  frequency  and  coupling  coefficient  and  measured  data 
supplied  by  the  manufacturer.  In  addition,  equations  for  the 
electrical  admittance  for  this  transducer  in  air  were 
developed  based  on  an  analytical  approach  published  by  Hong- 
zhang  Wang  (J.  Acoust.  Soc.  Am.  79,  164-176  (1986)).  A  Matlab 
program  was  written  to  implement  these  equations  and  to 
calculate  various  electromechanical  network  parameters.  The 
calculated  values  agreed  with  measured  values  to  within  5 
percent . 

An  in- water  harmonic  analysis,  in  which  the  transmitting 
voltage  response  and  the  directivity  pattern  are  computed,  was 
performed  using  the  ATILA  code.  In  general,  results  obtained 
using  the  finite  element  model  and  manufacturer's  measured 
values  agree  within  10  percent.  The  effects  of  changes  in 
material  properties  on  acoustic  performance  of  the  transducer 
were  investigated;  the  results  are  discussed. 
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I.  INTRODUCTION 


Free- flooded  piezoelectric  ring  transducers  leave  found 
increasing  use  in  deep  ocean  sonar  applications.  Though  simple 
in  concept,  the  ability  to  reliably  predict  the  performance  of 
such  devices  is  complicated  by  the  nature  of  their 
construction.  In  particular,  the  acoustic  performance  of  such 
transducers  depends  upon  the  physical  properties  of  the 
encapsulant  materials  used,  making  their  performance  difficult 
to  model . 

Standard  multiport  network  analysis  techniques,  long  used 
to  model  the  electromechanical  behavior  of  piezoelectric 
transducers,  are  usually  restricted  to  one -dimensional  motion. 
Newer  types  of  transducer  structures,  such  as  the  free- flooded 
ring,  that  have  more  than  one  important  degree  of  freedom,  do 
not  lend  themselves  to  this  technique.  Consequently ,  finite 
element  analysis  (FEA)  techniques,  which  provide  the 
capability  to  analyze  engineering  problems  of  many  degrees  of 
freedom,  are  frequently  employed. 

The  objectives  of  this  research  were  (1)  to  numerically 
model  a  free- flooded  piezoelectric  ring  sonar  transducer  and 
(2)  using  the  developed  model,  predict  the  changes  in  acoustic 
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performance  due  co  changes  in  the  elastic  properties  of  the 
encapsuiants  used. 

The  electroacoustic  properties  of  this  free- flooded 
piezoelectric  ring  sonar  transducer  were  numerically  modeled 
using  the  finite-element  code  "ATILA",  developed  by  the 
Institut  Superieur  a' Electronique  du  Nord  (Lille,  France) 
[Ref.  1].  Since  the  ring  is  constructed  of  ceramic  staves 
cemented  together,  the  effective  material  properties  for  the 
ring  were  determined  from  in-air  modal  analyses  using  ATILA. 
The  material  property  values  were  adjusted  to  obtain  the 
closest  agreement  between  the  calculated  ring  assembly 
resonance  frequency  and  coupling  coef ficient  and  measured  data 
supplied  by  the  manufacturer. 

In  addition  to  the  finite  element  analysis,  equations  for 
the  electrical  admittance  of  the  thin  ring  transducer  when 
unloaded  were  developed  based  on  an  analytical  approach 
published  by  Hong- zhang  Wang  [Ref.  2]  .  A  MATLAB  program  was 
written  to  implement  and  compute  various  electromechanical 
network  parameters.  The  results  obtained  from  the  analytical 
method  were  compared  with  the  measured  values  and  those 
obtained  by  the  finite  element  code. 

The  finite  element  model  built  was  further  used  to  predict 
the  changes  in  acoustical  properties  due  to  a  change  in  the 
type  of  encapsulant  materials  used  in  manufacturing. 

The  remainder  of  this  thesis  is  divided  into  seven 
chapters.  Chapter  II  describes  the  transducer  under  study. 
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Chapter  HI  gives  a  general  description  of  the  finite 
analysis  and  the  standard  istworx  analysis  if  piezo- 


t_  *0.  ryio  n  ^ 


ring  transducers.  Chapter  IV  describes  tne  development  or  the 
2  -  dimensional  finite  element  model  of  the  transducer  tor  use 
with  the  ATILA  code.  Chapter  7  presents  and  discusses  -he 
results  of  the  in-air  modal  analyses  and  m- water  narrr.cn ic 
analyses  using  ATILA,  as  well  as  the  results  cf  tne  analytical 
model.  Chanter  VI  presents  the  conclusions.  Chapter  VII 
provides  suggestions  for  future  work.  Appendix  A  contains  a 
cooy  of  the  incut  data  file  for  the  finite  element  mesh  used 
in  the  in-water  harmonic  analyses.  The  MATLAB  program  written 
to  implement  the  analytical  model  for  the  transducer  in  air  is 
attached  as  Appendix  E.  Appendix  C  provides  a  list  of  the 
typical  values  of  material  constants  used  in  the  modeling. 
Appendix  D  is  a  sample  calculation  of  the  effective  elastic 
constants  of  a  homogeneous  ceramic  ring  incorporating  the 
fiberglass  wrap. 
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II.  TRANSDUCER  DESCRIPTION 


Figure  1  shows  a  schematic  iiagram  :£  i  free-  flooded 
piezoelectric  ring  sonar  transducer.  When  immersed  m  water 
and  vibrating  radially,  such  a  transducer  exhibits  two 
fundamental  resonant  modes  chat  can  radiate  acoustically  into 
the  water.  They  are  the  radial  resonance  of  the  ring  and  the 
cavity  resonance  of  the  enclosed  water  column.  Radial  motion 
of  the  cylinder  walls  will  excite  the  symmetrical  cavity  modes 
in  the  enclosed  water  column.  The  radiai  resonance  is  the  more 
important  resonance. 


Figure  1:  The  Free  Flooded  Ring 
Transducer 


2  MODES:  CAVITY  &  RING 


Figure  3 :  The  Ring  transducer 


The  ceramic  plates  are  poled  tangentially  and  electrically 
wired  in  parallel.  Thus,  it  is  apparent  that,  during 
operation,  the  electric  field  will  either  simultaneouly  be 
parallel  or  antiparallel  to  the  poling  direction  shown  in 
Figure  3.  The  ring  is  wrapped  on  the  outside  by  epoxy- 
impregnated  fiberglass  of  0.03  in  (0.762  mm)  thickness,  which 
provides  a  comp-essive  stress  bias. 

A  cutaway  view  is  as  shown  in  Figure  4 . 
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Figure  4:  Cutaway  View  of  Ring  Transducer 

The  piezoelectric  ring  and  fiberglass  wrap  is  enclosed  by 
two  encapsulant  materials:  the  inner  one  is  a  polyurethane  and 
the  other  is  a  neoprene.  Of  particular  importance  is  the 
annular  mounting  plate  at  the  bottom,  which  has  an  inner 
radius  less  than  half  the  radius  of  the  transducer.  Because  of 
the  narrow  size  of  the  center  hole  of  the  mounting  plate,  it 
is  expected  that  the  mounting  plate  will  influence  the 
acoustic  cavity  mode  of  the  enclosed  water  column  in  the 
middle  of  the  transducer. 


7 


Appendix  C  lists  the  physical  dimensions  and  constants  for 
the  materials  of  the  transducer. 

The  main  operational  characteristics  of  the  transducer  as 
provided  by  the  manufacturer  are: 

a.  Resonance  frequencies  are  1.0  kHz  \ cavity)  and  2.5 
kHz  ring) ; 

b.  Source  Pressure  Level,  which  corresponds  to  the 
effective  pressure  on  the  acoustic  axis,  at  ring  resonance 
is  132  dB  ref  1  /iPa/V  at  1  meter; 

c.  Operational  Depth,  which  corresponds  to  the  maximum 
depth  where  the  performance  of  the  transducer  is  not 
compromised  -  unlimited. 

The  measured  admittance  circles  of  the  ceramic  ring 
assembly  without  and  with  fiberglass  wrapping  in  air  are  as 
shown  in  Figures  5  and  6,  respectively. 
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Figure  5:  Measured  Admittance  Circle  of  the 
Unwrapped  Ceramic  Ring  Assembly,  In 
Air 
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Figure  6:  Measured  Admittance  Circle  of  the 
Ring  Assembly  with  Fiberglass 
Wrap,  In  Air 
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The  measured  transmitting  voltage  response  TVR)  with  1 
V  drive  and  the  vertical  directivity  pattern  at  2  kHz  are  as 
shown  in  Figures  7  and  6,  respectively. 
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III.  METHODS  OF  ANALYSIS 


A.  INTRODUCTION 

This  chapter  presents  the  simple  approximate  expressions 
for  the  ring  and  cavity  mode  frequencies  from  G.W.  McMahon 
(Ref.  3]  and  discusses  the  theoretical  background  behind  the 
finite  element  code,  ATILA,  as  well  as  the  analytical 
approach,  suggested  by  Hong- zhang  Wang  [Ref.  2]  ,  to  obtain  the 
approximate  electrical  admittance  equations  for  the 
tangentially  polarized  piezoelectric  thin  circular  cylindrical 
tube  in  the  unloaded  case. 


B.  RING  AND  CAVITY  MODE  FREQUENCIES  [Ref.  3] 

The  behavior  of  the  cavity  and  ring  resonance  in  small 
length- to- radius  ratio  open  tubes  of  lead  zirconate  titanate 
ceramic  immersed  in  water  was  studied  by  G.W.  McMahon. 

When  the  tube  is  short  i.e.,  h/a  <§  tt  ,  G.W.  McMahon  gives 
the  following  formulas : 


(1) 
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where  u.  is  the  ring  mode  angular  frequency,  w  is  the  angular 
frequency  of  the  first  cavity  mode,  a  is  the  inside  radius,  h 
is  the  height  of  the  ring,  3  is  the  bulk  modulus  of  water,  t 
is  the  wall  thickness  of  ring  and  c  is  the  speed  of  sound  in 
water.  c0  is  the  velocity  of  sound  in  the  water  column  given 
by  Equation  5  below  and  a,  an  end  correction,  is  given  by  the 
empirical  formula: 

cc  =  0.633  -  0.106  Q  (3) 


The  dimensionless  frequency  parameter,  Q  -  co.a/c,,,  is  given  by 

Q(JL.  0 . 633  |  -0.106C2  =  1 
\  2  a  /  2 


The  value  of  Q  is  between  0.33  and  3.3.  The  velocity  of  sound 
within  the  tube,  c0,  is  effectively  less  than  the  velocity  of 
sound  c  in  open  water  because  of  the  finite  stiffness  of  the 
tube  wall .  It  is  given  by 


c0  =  c 


!l  + 


2JBa 

Y33t 


(5) 


The  calculated  value  of  the  cavity  mode  frequency  for  the 
transducer  is  presented  in  Chapter  V,  Section  C. 
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C.  FINITE  ELEMENT  ANALYSIS,  THE  ATILA  CODE 

The  application  of  finite  element  analysis  FEA;  to  solve 
boundary  value  problems  consists  of  the  transformation  of  the 
governing  differential  or  integral  equation  is;  into  a  multi - 
nodal  matrix  equation,  the  solution  of  which  represents  the 
discretized  solution  of  the  problem.  There  are  many  techniques 
to  obtain  a  finite  element  formulation  [Ref.  5,6, 7, a). 

ATILA  is  a  finite  element  code  developed  at  i'Institut 
Superieur  d' Electronique  du  Nord  (ISEN)  in  France  for  the 
analysis  of  underwater  transducers.  It  utilizes  the 
variational  formulation  of  the  finite  element  problem  [Ref. 
9,10,11,12,13]  . 

ATILA  uses  quadratic  isoparametric  elements.  Isoparametric 
means  the  same  polynomial  (quadratic  in  this  case)  is  used  to 
interpolate  both  the  geometry  and  the  field  variation.  ATILA 
has  46  different  types  of  elements,  including  shell,  plate, 
transition,  spring,  trilaminar,  and  two-  and  three-dimensional 
solid  elements  of  various  types.  With  ATILA,  it  is  possible  to 
model  lossy  elastic,  piezoelectric,  magnetostrictive,  magnetic 
and  composite  materials,  fluid,  solid- fluid  interfaces, 
dampers,  and  radiation  conditions. 

ATILA  can  perform:  (l)  static  analyses,  (2)  modal 
analyses,  which  correspond  to  a  free  vibration  problem,  where 
the  eigenfrequencies  and  eigenmodes  are  computed,  and  (3) 
harmonic  analyses  of  radiation  or  scattering  problems,  which 
correspond  to  a  forced  vibration  problem,  the  excitation  being 
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the  voltage  applied  across  the  electrical  terminals  of  cne 
transducer  or  external  forces  applied  to  the  nodes. 

1.  Harmonic  Analysis  of  a  Radiating  Piezoelectric 

Transducer 

This  problem  is  governed  by  the  equation  of  motion  in 
the  elastic  and  piezoelectric  structures,  by  Poisson's 
equation  in  the  piezoelectric  structures,  and  by  Helmholtz's 
equation  in  the  fluid.  Appropriate  boundary  conditions  are 
defined,  both  on  the  solid- fluid  interface  and  on  the  external 
fluid  boundary,  which  must  simulate  the  appropriate  radiation 
condition. 

The  solid  equation  of  motion  is  given  by  [Ref.  12] : 


P 


dzUj 

W 


a  a, 

dx, 


(6) 


where  p  is  the  solid  material  density,  u  is  the  displacement 
vector,  t  is  time,  a  is  the  stress  tensor,  and  x,-  is  a 
coordinate  direction.  Here  i  and  j  can  be  1,  2  and  3,  and  the 
Einstein  notation  is  used,  where  summation  is  implied  over 
repeated  indices  in  the  same  term. 

Poisson's  equation  is  given  by  [Ref.  12]: 


dDt 

Tx, 


=  0 


(7) 
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where  D  is  che  electric  displacement  vector  and  is  a 
coordinate  direction;  i  can  be  1 .  2  and  3. 

The  linearized, lossless  Helmholtz  Equation  for  the 
propagation  of  sound  in  fluids  is  given  by  [Ref.  14] : 

V2  p  -  EE  =  0  (8) 

c~  at; 

where  V2  is  the  three-dimensional  Laplacian  operator,  p  is  the 
acoustic  pressure,  and  t  is  time. 

In  piezoelectric  materials  the  stress  tensor  and  the 
electric  displacement  vector  can  be  related  to  the  strain 
tensor  and  electric  field  vector  via  material  properties  by 
the  following  constitutive  equations,  which  neglect  magnetic 
and  pyroelectric  effects  [Ref.  12]: 

%]  =  c* >ju  ~  eyk  ^  ^ 

D,  .  ea  S„  .  Is,  Bj  (10) 

where  (21  is  the  stress  tensor,  [S]  is  the  strain  tensor,  E  is 
the  electric  field  vector,  D  is  the  electric  displacement 
vector,  [c®]  is  the  constant  electric  field  elastic  stiffness 
tensor,  [e]  is  the  piezoelectric  tensor,  and  [e5]  is  the 
constant  strain  dielectric  tensor;  i,  j,  k  and  1  can  be  equal 
to  l,  2  and  3. 
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Ultimately  the  solution  must  be  expressed  in  terms  c: 
displacements  and  electric  potentials.  To  this  end  the 
following  two  equations  from  elasticity  ana  electricity, 
respectively,  are  used  'Rer.  12]: 


S, 


1  (  du  du 

2  '  dx  Ox. 


(11) 


B, 


d<$ 

Jx: 


(12) 


where  [S]  is  the  strain  tensor,  u  is  the  displacement  vector, 
Xj  is  a  coordinate  direction,  E  is  the  electric  field  vector 
and  <P  is  the  electrical  potential;  i  and  j  can  be  equal  to  1, 
2  and  3.  It  must  be  noted  that  the  quasi-static  approximation 
[Ref.  4}  is  used  in  writing  Equation  12  for  the  low  (acoustic) 
frequencies  involved  here. 

The  boundary  conditions  and  prescribed  excitations  at 
each  node  can  be  defined  either  by  a  displacement  or  an 
applied  force,  an  electrical  potential  or  an  electrical 
charge,  or  an  acoustic  pressure. 

In  ATILA,  the  previous  seven  equations  are  transformed 
into  the  following  matrix  equation  [Ref.  1]  : 

[FJ-ah[M]  [K^]  ~[L) 

[Kui]r  [**]  [0] 

-p2c2u2[L]r  [0]r  [H]  ~ur  [Mj  ] 


u 

F 

* 

-<7 

p 

p2c2\j/ 

- 

L.  -J 
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wnere : 


U:  vector  of  the  nodal  values  of  the  components  of 

the  displacement  field, 


vector 

of  the  nodal  values 

the  electrtai 

potenti 

al , 

P: 

vector 

of  the  nodal  values  of 

the 

treasure  field, 

F: 

vector 

of  the  nodal  values 

of 

the  externally 

applied 

forces , 

g: 

vector 

of  the  nodal  values 

of 

the  externally 

applied 

electrical  charges, 

vector 

of  the  nodal  values 

of 

the  integrated 

normal 

derivative  of  the 

externally  applied 

pressure  on  the  surface  boundary  S, 

[JCUU]  :  Finite  Element  stiffness  matrix, 

[Jfu4)j  :  Finite  Element  pieozoelectric  matrix, 

[Kqq]  •  Finite  Element  dielectric  matrix, 

[M]  :  Finite  Element  consistent  mass  matrix, 

[H]  :  finite  Element  fluid  pseudo-  stiffness  matrix, 

[M,]  :  Finite  Element  consistent  fluid  pseudo-  mass 

matrix, 

[L] :  Finite  Element  coupling  matrix  at  the  fluid 

structure  interface, 

[0] :  zero  matrix, 

«:  angular  frequency, 

p:  fluid  density, 
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c:  fluid  sound  speed. 

*:  means  transposed. 

The  results  of  this  analysis  for  each  input  frequency 
are  the  complex  displacement,  rotation,  and  electrical 
potential  fields  at  each  transducer  node,  the  complex  pressure 
field  at  each  fluid  node,  and  the  complex  electrical  impedance 
and  admittance. 

2.  Modal  Analysis  of  a  Piezoelectric  Transducer 

This  problem  is  governed  by  the  equations  or  motion  m 
the  elastic  and  piezoelectric  structures,  and  by  Poisson's 
equation  in  the  piezoelectric  structures.  The  matrix  equation 
governing  this  problem  is  easily  obtained  from  that  aescnbed 
in  the  previous  section.  In  a  modal  analysis  there  is  no  fluid 
and  there  are  no  external  forces  applied  (the  natural  boundary 
conditions) ,  so  the  third  row  and  column  of  Equation  13  become 
irrelevant,  and  F  is  replaced  by  0,  resulting  in 


[KJ  -ur[M] 

~U~ 

"  0 

$ 

where  the  elements  are  defined  by  Equation  13. 

In  this  equation  the  resonance  condition,  which 
corresponds  to  the  electrical  short-circuit  condition,  is 
obtained  by  setting  $=0 .  The  anti -resonance  condition,  which 
corresponds  to  the  electrical  open- circuit  condition,  is 
obtained  by  setting  g=0  [Ref  15] . 
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The  results  of  this  analysis  are  the  eigenf requencies 
and  eigenmoaes.  The  maximum  number  of  modes,  which  must  be 
specified  by  "he  user,  is  IOC  . 

D.  STANDARD  NETWORK  ANALYSIS  [Ref.  2] 

The  state  equations  for  piezoceramics  as  stated  m 
Equation  S  and  10  can  generally  be  expressed  as  follows: 

=  SeIJ  T.  -  d.  R  15! 

D<  =  d_.  T  -  eTu  E,  '16) 

where  [  T]  is  the  stress  tensor,  [S]  is  the  strain  tensor,  E 
is  the  electric  field  vector,  D  is  the  electric  displacement 
vector,  [a5]  is  the  constant  electric  field  elastic  compliance 
tensor,  [d]  is  the  transposed  piezoelectric  stress  tensor,  and 
[e  /  is  the  constant  stress  permittivity  tensor;  i,  j,  k  and 
1  can  be  equal  to  1,  2  and  3.  The  matrix  notation  has  been 
used  [Ref .  16]  . 

If  it  is  assumed  that  the  ceramic  is  isotropic  in  the 
plane  transverse  to  the  poling,  the  piezoelectric  equations  in 
the  cylindrical  coordinate  system  may  be  written  as 


20 


where  sn~s33,  s,v 

U  =  (Ur,  Ug,  Uz 


respectively . 


5  -  S  - 

rr  .  \  ■ 

c  -  cy£  /r' 


;  «*•  vfj  A  n  ¥ 

3Eh 

S.. = s  “ , :  Trr  -  s  ",  s  ^  ,  |  x ..  ~-o.X{Eg 


SrD~s  ^-d^E, 


rV 


*15* 


S.t=sEuT.t~d'5E. 
q  =o£  T 

°r.  °  66 ■‘■ft: 


Dr~d\ sTr1)~  €Tnr.r 

l^=c^3i  Trr-rdgjTgg-rd}i T„ -  e 
^:  =  di;T:tl~  €  t|-fi. 


117) 


=s13.  If  we  let  the  displacement  be  expressed  as 
> ,  then  the  strain  components  are 


S  = 


S«=* 


dur 

Ur 

duf_  u9  i  dur 

tT  ~W 


«~Tr 
1  du, 


S.„  = 


u. 


r  W*T 
du9  i  3u. 

5z  ~r  d6 
du. 


(18) 


_  Bu.  dur 

S*'UF  +  Ui" 
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The  equations  of  motion  of  the  elastic  medium  in  the 


cylindrical  coordinate  system  are 

.  dT„r  :  dT„  dT„  7,t  - 

pu''"Tr  "  r “9T 
..  dT^  :  dZH  dT„  2  T, 

pu*~  TF  "  r  UF  '  "<7F  *  T~ 

-  3T.  !  ar,.  3Tn  Tn 

pU:~-57~^-JF^-ZF~~T 


1.  Vibration  of  Unloaded  Tube 

For  the  tangentially  polarized  ring,  the  following 
assumptions  were  made: 

a.  the  electric  field,  E  is  approximated  as 

(0,  0} 

b.  the  unloaded  thin- wall  tube  vibrates  with 
axial  symmetry 

c.  the  stress  components,  Txz,  Tzr,  Trr,  T^, 
and  TQr  on  the  outer,  inner  surfaces  and  in 
the  whole  body  of  the  tube  are  neglected. 

d.  Tzz  =  o  at  z  =  ±  1/2 

e.  Ug  =  0; 

f.  dT^/dd  =  dTQe/de  =  dTza/dQ  =  0 


22 


With  the  assumptions  mentioned,  the  piezoeiec; 
equations  may  cs  simplified  as 


~  L  '  1  ''T'  '  .  {S/M  '  J^-_z  r  I 

T-.  =  [  Y£7  ( 1  -  j:;  (<r  (d„  -<Jd,3)  E„  | 

n  —  W  'T1  _  W  T*  —  *:  •*"  C 

—  ft  -*'33  ^  “  13  *-J« 


where  Y£  =  l/sMs  is  the  Young's  modulus  and  <r  =  -s,;*/s,,A  is  the 
Poisson  ratio. 

The  vibration  equations  of  the  thin-walled 
piezoceramic  cylindrical  tube  can  also  be  simplified  as 


pu.  =  -  — 

r 

pue  =  0 

.  dT„ 

pu:  =  -jj- 


By  combining  the  piezoelectric  equations  and  the  equation 
of  motion,  the  following  acoustic  wave  equation  may  be 


obtained: 


d2u.  , 

- :  +  Kju.  =  0 

dzz  ‘  ‘ 


where 


K2  =  — 
c„ 


2  fa(w) 
[(u/ur)2- 
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,u )l 


24) 


and  ojr  =  c./a.  c.  is  the  acoustic  wave  speed  an  the  ceramic 
given  by: 

(25) 


Applying  the  boundary  conditions,  we  obtain  the  following 

u.  =  Ag  sin (.KjZ)  (26) 

^  =  JaEt  (27) 


*4^  (U)+<r2]  cos 


If  a  and  b  are  the  inside  and  outside  radii  of  the  ring 
and  N  is  the  number  of  silvered  electrodes,  the  total  electric 
charge  on  the  electrodes  is  given  by 

Q  =  NQX  =  J*D#  drdz  (29) 
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The  electric  potential  between  the  two  near  electrodes  is 
%  and  from  Maxwell’s  equation,  it  follows  that 


Combining  the  aforementioned  relations,  the  admittance 
equation  are  obtained  as 


f.(oj) 


“ 

- 

(f ,  ( ui)  {l*aT)^a((J  +  T\y 

tan  (  —K,l) 

2  J 

!.  a  +  r ) : 

<1  -<r);  [  (to/ a?.;  * -1] 

±K,1 

1  -cr 

2  ' 

1 

J 

(31; 


-jtoC 


where 


V  = 


:  33 


C,  *  C0[l-(VV(l-<72)  >  <1+2<XT  +  T2)  3 
In  (jb/a) 


c0 


27T 

T  ~  (  ^3/  ) 


These  equations  are  applied  in  obtaining  the  results  in 
Chapter  V,  Section  E. 
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IV.  TWO-DIMENSIONAL  FINITE-ELEMENT  MODEL 

A.  INTRODUCTION 

In  order  to  efficiently  model  the  free  flooded  ring,  as 
shown  in  Figure  1,  the  structure  may  be  considered  to  be  an 
ax i symmetric  body.  Appropriate  modifications  to  the  material 
properties  must  be  made  to  model  the  structure  as 
axisymmetric .  A  continuing  problem  in  the  analysis  of  these 
piezoelectric  devices  is  the  determination  of  these  material 
properties . 

A  two-dimensional  finite  element  model  was  developed  using 
8 -node  isoparametric  quadrilateral  elements  for  the 
piezoelectric  ceramic,  the  elastic  encapsulant  materials,  the 
mounting  plate,  supporting  structures,  and  the  fluid.  The 
fluid  domain  extends  to  lm  from  the  acoustic  center.  6 -node 
isoparametric  linear  elements  were  used  for  the  solid- fluid 
interface  and  3 -node  linear  elements  were  used  for  the  dipolar 
radiating  damper.  These  were  used  to  terminate  the  fluid  mesh. 
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B.  CHARACTERISTICS  OP  THE  MODEL 
1.  Material  Properties 

a.  Unwrapped  Piezoelectric  ceramic  Ring 

As  described  before,  the  transducer  has  72 
tangentially  poled  lead  zirconate  titanate  ceramic  staves, 
glued  together  with  epoxy  resin,  and  arranged  in  a  0.3668  m 
diameter  ring. 

In  order  to  simplify  the  finite  element  model,  an 
equivalent  homogeneous  piezoelectric  ceramic  ring  with 
material  properties  determined  by  a  combination  of  the 
piezoelectric  ceramic  and  epoxy  resin  properties  was  used. 
Because  the  appropriate  combination  of  properties  is  not 
known,  "smeared"  material  properties  for  the  homogeneous  ring 
were  obtained  by  adjusting  the  circumferential  compliance  so 
that  the  resonance  frequencies  obtained  using  the  ATILA  code 
tolerably  matched  the  in- air  experimental  measurements  of  the 
resonance  frequencies  for  the  unwrapped  ring. 

A  limitation  of  an  axisymmetric  ATILA  model  is  that 
the  tangential  poling  of  the  segmented  ring  cannot  be  modeled 
directly.  It  has  to  be  modeled  as  a  radially-poled  ring  as 
illustrated  in  Figure  9 . 
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MODELED  AS  RADIALLY-POLED  RING 


Figure  9:  Modeled  as  a  Radially- Poled  Ring 

The  transformation  from  the  tangentially  to  radial 
polarization  is  accomplished  by  a  suitable  exchange  of  the 
elastic,  piezoelectric,  and  dielectric  tensors  as  described  in 
the  ATILA  manual  [Ref.  1]  .  Since  the  potential  degree  of 
freedom  is  a  scalar,  the  change  in  polarization  direction  is 
simply  a  geometrical  transformation.  Thus,  to  model  the 
applied  potential  of  the  tangentially  poled  ring  using  the 
radially  poled  ring,  the  radial  potential  has  to  be  corrected 
by  a  factor  equal  to  the  ratio  of  the  radial  thickness  to  the 
stave  thickness,  d/e  in  Figure  9.  Similarly,  the  results 
computed  in  a  harmonic  analysis  for  the  radially-poled 
axisymmetric  model  have  to  be  modified  as  follows:  (1)  the 
displacement  and  pressure  magnitudes  must  be  divided  by  the 
ratio  of  the  circumferential  length  of  a  "smeared" 
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piezoelectric  element,  which  is  equal  to  1/72  of  the  ring 
circumferential  length,  to  its  thickness,  and  (2)  multiply  the 
electrical  impedance  by  the  square  of  the  same  ratio.  This 
ratio  is  2.602  for  the  transducer  under  study. 

Jb.  Wrapped  Piezoelectric  Ceramic  Ring 

Normally,  the  fiberglass  belt  wrapped  around  the 
ceramic  ring  (Figure  3)  is  modeled  as  thin  shell  element. 
However,  to  simplify  the  model  and  to  reduce  the  number  of 
degrees  of  freedom,  the  0.762  mm  thick  fiberglass  belt  in  our 
case  was  not  modeled  as  a  separate  element.  Instead,  the 
material  properties  of  the  equivalent  homogeneous  ceramic  ring 
were  modified.  A  sample  calculation  of  the  approximate  new 
material  constants  is  given  in  Appendix  D.  These  approximate 
values  were  then  used  for  a  modal  analysis  using  ATILA  and 
were  fine  tuned  by  adjusting  first  the  circumferential 
compliance  so  that  the  resonance  frequency  obtained  tolerably 
matched  the  measured  values  provided  by  the  manufacturer.  Then 
the  other  compliances  were  adjusted  to  obtain  a  close  match  of 
the  coupling  coefficient  to  the  measured  value.  This  scheme  is 
considered  sufficient  for  the  model  even  though  a  scheme 
implemented  by  McMahon  and  Armstrong  [Ref.  17]  ,  which  fits  the 
resonance  and  antiresonance  frequencies  and  the  electrical 
capacitance  taking  into  consideration  the  compression  bias 
effect  of  the  fiberglass  wrapping,  might  give  a  more  accurate 
result . 
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c.  Encapsulant  Materials  and  Mounting  Plate 

The  encapsulant  materials,  namely  the  polyurethane 
and  the  neoprene,  the  reflector  ring  thigh  carbon  steer) ,  and 
the  mounting  plate,  a  fiberglass -epoxy  composite,  were  modeled 
individually  as  8 -node  quadrilateral  elastic  elements.  The 
material  properties  were  provided  by  the  manufacturer. 

2 .  Types  of  Elements 

The  quadratic  isoparametric  elements  listed  m  Table 
1,  which  are  described  in  the  ATILA  User's  Manual  [Ref.  13, 
were  used: 

Table  1:  TYPES  OF  ELEMENTS  USED 


Region 

Element 

Geometry 

Piezoelectric  ring 

QUAD08P 

8 -node  quadrilateral 

Encapsulant  ( Boot  a , 
Spacer,  Reflector 
Ring)  &  Mounting 
Plate 

QUAD08E 

8 -node  quadrilateral 

Interface  solid- fluid 

LINE06I 

6 -node  isoparametric 

Fluid 

QUAD08F 

8 -node  quadrilateral 

Radiation  surface 

LINE03R 

3 -node  linear 

30 


3.  Constraints  on  Mesh  Design  [Ref.  1] 

Design  of  the  mesh  was  guided  by  the  following 
constraints : 

a.  Aspect  ratio 

The  aspect  ratio  of  each  element  should  not  be 
greater  than  3,  although  4  is  considered  an  acceptable,  though 
less  conservative,  value. 

b.  Internal  angles 

The  angles  between  adjacent  sides  of  quadrilaterals 
should  be  between  45  and  135  degrees,  although  30  degrees  and 
150  degrees  are  considered,  respectively,  acceptable,  though 
less  conservative,  values. 

c.  Element  size 

As  ATILA  utilizes  quadratic  interpolation 
functions,  the  size  of  each  element  must  not  be  greater  than 
one -fourth  wavelength  at  the  highest  frequency  of  interest. 

d.  Interelement  compatibility 

The  mesh  should  be  built  in  such  a  way  that 
adjacent  elements  have  adjoining  sides  with  collocated  nodes 
to  ensure  accurate  interpolation  at  their  interfaces. 

e.  Radiation  boundary  elements 

For  in-water  harmonic  analyses  (radiation  problems) 
the  fluid  mesh  outer  limit  must  be  spherical.  This  is  required 
by  the  radiation  damping  elements  available  in  the  ATILA  code . 
ATILA  offers  damper  elements  that  are  attached  to  the  external 
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surface  of  the  surrounding  fluid  domain  and  are  designed  :o 
absorb  various  components  of  the  radiated  multipolar  acoustic 
field  to  overcome  the  problem  of  modeling  an  infinite  fluid 
domain.  The  ATILA  offers  the  monopoly  and  dipole  radiation 
damping  elements.  The  latter  includes  not  only  the  monopole 
term  of  the  radiated  field  multipolar  expansion,  but  also  the 
dipole  term.  Dipolar  damping  elements  were  selected  to 
terminate  the  fluid  mesh  because,  in  addition  to  providing  a 
more  accurate  solution,  a  smaller  mesh  can  be  employed, 
resulting  in  a  lower  computational  cost. 

A  fluid  mesh  outer  limit  radius  greater  than  the 
far- field  distance  is  desirable  to  compute  the  acoustic  source 
pressure  level  and  to  compare  compu’  ed  and  measured  acoustic 
pressure  data.  The  far  field  crite:  ia  for  a  baffled  rigid 
piston- like  source  is  given  by  [Ref.  15} : 


and 


Rff  >  D 


(34) 


where  D  is  the  diameter  of  the  ring  which  is  0.3795  m. 

The  minimum  wavelength,  X,  corresponding  to  the 
maximum  frequency  of  interest  is  given  by: 


X  =  ~  =  0.375m 


(35) 
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where  is  Che  speed  of  sound  in  water,  1500  m/s,  and  f  is  the 
highest  frequency  of  interest.  4000  Hz,  giving  Rff  >  0.3  84  m. 

The  boundary  was  placed  at  a  radius  R  equal  to  1  m 
from  the  transducer's  acoustical  center,  which  is  more  than 
2.5  times  the  far -field  limit  of  the  equivalent  piston- like 
source  at  the  resonance  frequency  [Ref.  18} . 
f.  Boundary  Conditions 

Boundary  condition  data  entries  enable  the  user  to 
force,  for  example,  clamped,  hinged,  or  simply- supported 
conditions,  or  to  set  the  master  degrees  -  of  -  freedom.  These 
conditions  are  generally  associated  with  symmetry  planes  or 
axes,  electroded  or  pressure  release  surfaces. 

In  a  harmonic  analysis  of  a  piezoelectric 
structure,  the  degree  -  of  -  freedom,  D,  for  all  the  nodes 
associated  with  the  applied  potential,  V,  to  the  plane  or 
line,  P,  containing  the  node  N  has  to  be  a  master  degree  of 
freedom  and  identical.  This  is  realised  by  using  the  following 
line  in  our  data  file: 

-N  -D  -P  V 

for  example, 

-1  -4  -2  1 

where  the  degree  of  freedom  4  indicates  electrical  potential 
on  a  plane  2  (normal  to  the  OY  axis)  passing  through  the  node 
number  1  with  an  excitation  voltage  of  l  V.  The  negative  sign 
in  front  of  the  2  indicates  that  the  electrical  potential 
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degree  of  freedom  of  ail  nodes  on  chat  plane  are  identical, 
i . e . ,  IV. 

Since  our  model  is  a  2-D  model,  the  degree -of  - 
freedom  in  the  plane  normal  no  the  OZ  axis  are  ail  deleted 
using  the  following  data  line: 

-N  D  ? 

where  D  is  3  which  corresponds  to  the  displacement  in  the  z- 
direction  and  F  is  3  which  corresponds  to  the  plane  normal  to 
the  OZ  axis. 

4.  Final  Mesh  Designs 

With  the  above  considerations,  an  in- air 
axisymmetricai  mesh  about  the  OX  axis  was  designed  for  fitting 
the  material  properties.  Only  half  of  the  ring  cross-section 
needs  to  be  modeled.  The  bare  ceramic  staves  without  the 
fiberglass  wrap  were  divided  into  five  elements  to  meet  the 
aspect  ratio  requirement  of  the  ATILA  code.  The  in-air  model 
is  as  shown  in  Figure  10. 
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IN-AIR 
F.  EE .  MODEL. 


x 


Figure  10:  In- Air  Model 


Using  the  pre-processor  mesh  generator  MOSAIQUE  [Ref.  1] , 
the  in-water  mesh  without  the  mounting  plate  as  shown  in 
Figure  11  was  generated.  The  transducer  mesh  consists  of  40 
solid  elements,  21  interface  elements,  166  fluid  elements,  ant 
15  radiating  elements,  totaling  2884  DOF.  This  mesh  was  also 
used  to  investigate  effects  of  encapsulant  material  changes. 

The  final  mesh  to  be  developed  was  the  in-water  model  with 
the  mounting  plate.  The  total  mesh  is  shown  in  Figure  12  and 
an  enlarged  view  of  the  solid  mesh  with  the  mounting  plate  is 
as  shown  in  Figure  13.  There  were  a  total  of  87  solid 
elements,  48  interface  elements,  406  fluid  elements  and  35 
radiating  elements  with  a  total  of  6600  DOF.  The  input  data 
file  of  this  final  mesh  is  given  in  Appendix  A. 
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Figure  11:  In- Water  F.E.  Model  Without  Mounting  Plate 
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IN-WATER 
F.E.  MODEL. 
WITH 

MOUNTING 

PLATE 


Figure  12:  In-Water  Model  With  Mounting  Plate 
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Figure  13 :  The  Enlarged  View 
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V.  RESULTS 


A.  MATERIAL  PROPERTIES  OF  EQUIVALENT  HOMOGENEOUS  RING 

The  material  properties  of  the  equivalent  homogeneous 
piezoelectric  ceramic  ring,  unwrapped  and  wrapped,  were 
determined  by  modal  analysis  using  ATILA.  This  analysis 
corresponds  to  a  free  vibration  problem  in  which  the 
eigenf requencies  and  eigenmodes  are  computed. 

1.  Unwrapped  Equivalent  Homogeneous  Ring 

First,  a  simple  transducer  model  of  only  the  bare 
ceramic  ring  without  the  fiberglass  wrap,  as  shown  in  Figure 
10,  was  designed.  A  modal  analysis  of  this  model  was  conducted 
using  book  values  for  the  piezoelectric  properties.  Table  2 
lists  the  calculated  resonance  and  antiresonance  frequencies 
and  the  coupling  coefficient  for  the  first  ten  modes  of 
vibration.  The  coupling  coefficient  is  calculated  using  the 
relationship  [Ref.  19]: 


k2 


(36) 


where  k  is  the  coupling  coefficient,  fr  is  the  resonance 
frequency  and  fa  is  the  antiresonance  frequency. 
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TABLE  2:  RESULTS  OF  THE  MODAL  ANALYSIS  OF  UNWRAPPED 


PIEZOELECTRIC  RING  USING  SOOK  VALUES  FOR 
THE  MATERIAL  PROPERTIES 


No. 

RESONANCE 

ANTIRESONANCE 

COUPLING 

(  HZ  ) 

(  HZ  } 

COEFF . ,k 

1 

o.  . 

1 . 0  60 ” j 

0 . 86e-  3 

-  0 

2  . 

2491 

3287 

0.652 

3. 

2499 

2499 

0 

4. 

— — - - 

2611 

2611 

0 

5. 

3287 

3488 

0.335 

6. 

4889 

4889 

0 

7. 

7393 

7393 

0 

8. 

10402 

10446 

0.092 

9  . 

10728 

10728 

0 

10  . 

14900 

14900 

0 

The  ring  mode  of  vibration  (mode  2)  is  the  most 
important  in  operation  and  is  shown  in  Figure  14 .  The  dashed 
lines  in  Figure  14  correspond  to  the  rest  position,  and  the 
deformed  shape  is  in  solid  lines. 


To  obtain  the  smeared  properties  of  the  unwrapped  ring 
assembly,  the  circumferential  compliance  s33  was  adjusted  so 
that  the  ring  mode  resonance  frequency  obtained  for  the 
unwrapped  ceramic  ring  in  air  using  the  ATILA  model  matched 
the  manufacturer's  measured  values. 

Table  3  compares  the  results  obtained  using  typical 
elastic  constant  values  provided  in  handbooks  and  the  adjusted 
values . 
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TABLE  3:  RESULTS  FOR  RING  MODE  OF  UNWRAPPED  CERAMIC  RING 


USING  BOOK  VALUES  AND  ADJUSTED  VALUES 


■ 

MEASURED 

USING  BOOK 

VALUES 

ADJUSTED 

VALUES 

Resonance 

2356  Hz 

2491  Hz 

2356  Hz 

(5.7%) 

f 0.0%) 

Anti resonance 

4224  Hz 

3287  Hz 

3177  Hz 

Coupling 

0.665 

0.652 

0 . 671 

Coefficient,  k 

(0.9%) 

Note:  The  number  in  brackets  are  the  percentage  error. 

A  comparison  of  the  value  for  the  ring  mode  frequency 
for  the  ATILA  2-D  in-air  model  for  the  unwrapped  ceramic  using 
the  adjusted  values  for  elastic  constants  and  the  in-air 
experimental  results  provided  by  the  manufacturer  shows  the 
expected  match  of  resonance  frequency.  The  computed  coupling 
coefficient,  however,  was  0.9%  higher  than  the  measured  value. 
The  difference  is  due  to  the  fact  that  only  the  elastic 
constant  s33  was  adjusted.  (Several  ATILA  modal  analysis  runs 
were  made  with  adjustments  made  to  the  elastic  constants  sM 
and  s12  just  to  see  the  effect  on  the  ring  resonance  frequency 
and  coupling  coefficient.  The  resonance  frequency  changed  only 


very  slightly  and  the  change  to  the  coupling  coefficient  was 
also  very  small.;  It  was  fait  that  a  0.9%  difference  is 
acceptable  for  our  purposes.  For  our  results,  the  value  of  sI3 
was  increased  by  12  percent  i.e.,  more  compliance;  from  the 
book  value  of  G.1539E-1Q  to  G.1724E-1G.  This  is  consistent 
with  our  expectation  that  a  ring  made  of  glued  segments  is 
more  compliant  chan  a  solid  piezoelectric  ceramic  ring.  This 
value  of  s33  was  used  for  calculating  the  first -guess  smeared 
material  properties  for  the  fiberglass  wrapped  ceramic  ring  as 
described  below. 

2.  Wrapped  Equivalent  Homogeneous  Ring 

Smeared  material  properties  for  the  homogeneous 
equivalent  wrapped  ceramic  ring  assembly  were  obtained  by  a 
four- step  process.  First,  the  coarse  values  for  the  elastic 
constants  were  calculated.  A  sample  calculation  is  given  in 
Appendix  D.  The  circumferential  compliance  s33  is  then  adjusted 
so  that  the  calculated  ring  resonance  frequency  matched  the 
measured  frequency.  The  other  elastic  constants  were  then 
adjusted  by  a  common  factor  such  that  the  computed  coupling 
coefficient  matched  that  of  the  measured  value.  The  results 
for  the  ring  resonance  frequency  and  coupling  coefficient  are 
given  in  Table  4: 
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TABLE  4:  RESULTS  OF  THE  WRAPPED  EQUIVALENT  CERAMIC  RING 


MEASURED 

USING  VALUES 

OBTAINED 

EARLIER 

ADJUSTED 

VALUES 

RESONANCE 

2462  Hz 

2442  Hz 

2465  Hz 

FREQUENCY 

(0.8%) 

(0.1%) 

COUPLING 

Q  .  646 

0.65 

0.644 

COEFFICIENT 

(0.62%) 

(0.3%) 

Note:  The  number  in  brackets  are  the  percentage  error. 

Table  5  lists  the  values  of  the  elastic  constants 
calculated  as  described  in  Appendix  D,  and  the  final  adjusted 
values.  These  latter  values  were  then  used  for  in-water 
harmonic  analyses: 
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TABLE  5:  ELASTIC  CONSTANTS  USED  FOR  IN -WATER  HARMONIC 


ANALYSES 


COMPLIANCE 

COMPUTED  AS  IN 

ADJUSTED  m*/N)  | 

l 

CONSTANTS 

APPENDIX  D  ,m;/N) 

1 

1 

Sjj®  S02 

1 . 142e- 11 

1 .3426-11 

Sn= 

-4 . 9e- 12 

-4 . OSe-12 

s13=s31 

-  5 . 3e- 12 

.J 

-4 . 43e- 12 

S03  = 

-  5 . 3e- 12 

-  4 . 43e- 12 

S33 

1 . 722e- 11 

1 . 702e- 11 

SiA~  S55 

3.9e-ll 

3 . 9e- 11 

S66 

3 . 25e- 11 

3 . 58e- 11 

The  other  material  constants,  such  as  the  dielectric 
constants  and  the  piezoelectric  constants,  were  not  adjusted. 
The  effective  density  of  the  homogeneous  ring  incorporating 
the  fiberglass  wrap  was  obtained  by  dividing  the  total  mass  of 
the  ceramic  and  fiberglass  by  the  volume  of  the  homogeneous 
ceramic  ring. 
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B.  IN-WATER  HARMONIC  ANALYSES 

A  harmonic  analysis  corresponds  to  a  forced  vibration 
problem,  the  excitation  being  the  voltage  applied  across  the 
electrical  terminals  of  the  transducer.  The  applied  voltage  in 
our  case  is  1  Vrms .  Internal  material  losses  for  the 
encapsulant  materials  were  included  in  this  model  {but  not  for 
the  ceramic  ring  assembly')  .  These  were  obtained  from  the  loss 
tangents  for  the  respective  materials.  Values  of  0.4  and  0.7 
were  used  for  the  loss  tangent  for  polyurethane  and  neoprene, 
respectively  [Ref.  20]  . 

1.  Mounting  Plate  Not  Included 

Using  the  material  property  values  obtained  from  in¬ 
air  modal  analyses,  described  earlier,  an  in- water  model 
without  the  mounting  plate  was  first  nalyzed.  The  mesh  is 
shown  in  Figure  11. 

The  plot  in  Figure  15  depicts  the  transmitting  voltage 
response  curve  obtained  by  ATILA  for  this  model  along  with  the 
corresponding  manufacturer's  data  for  the  actual  transducer 
(with  the  mounting  plate) . 

The  computed  transmitting  voltage  response  was 
generally  of  the  same  form  except  for  the  lower  frequency 
region  between  1000  to  1750  Hz.  The  discrepancy  in  this  region 
is  probably  due  to  not  including  the  mounting  plate  in  this 
model.  The  mounting  plate  is  expected  to  change  the 
transmitting  voltage  response  at  the  cavity  resonance  as  well 


dB  r< 


as  the  cavity  resonance  frequency  inseif .  The  TVR  ac  the  ring 
resonance  frequency  was  caicuiaced  to  be  aoout  2  a5  lower  than 
that  measured  at  the  same  rreauencv. 


Figure  15:  Comparison  Of  The  Measured  Transmitting  Voltage 
Response  With  That  Computed  For  In- water  Model 
without  Mounting  Plate 
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2.  Mounting  Plate  Included 

An  in- water  ATILA  model  including  the  mounting  plate 
was  then  developed.  The  mesn  are  shown  in  Figures  11  and  13. 
A  comparison  of  the  calculated  transmitting  voltage  response 
with  the  measured  response  is  shown  in  Figure  16  below: 


140  r 

i 


FREQUENCY 


Figure  16:  Comparison  Of  Measured  TVR  With  That  Calculated 
For  In-water  Model  with  Mounting  Plate 


The  agreement:  between  the  calculated  and  measured  TVR 
is  about  the  same  as  the  previous  model.  It  is  jreatly 
improved  m  the  cavity  mode  region  and  the  ring  resonance 
t  recxuency  calculated  is  -xactly  tne  same  as  the  measured 
values.  However,  there  is  a  deep  dip  at  150G  Hz.  If  this  dip 
can  be  made  shallower,  then  the  overall  agreement  would  be 
excellent.  It  is  hoped  that  the  results  will  be  better  when 
accurate  material  property  values  are  available  for  the 
polyurethane  and  neoprene  encapsuiants .  The  variation  in  TVR 
with  encapsulant  material  properties  is  explored  in  the 
Section  D. 

C.  COMPARISON  OF  CAVITY  MODE  FREQUENCY  [Ref.  3) 

The  results  of  the  previous  section  for  the  cavity  mode 
frequency  can  be  compared  with  the  approximate  analytical 
formulae  given  by  G.W.  McMahon  and  present  previously  in 
Chapter  III. 

Using  relationships  given  by  G.W.  McMahon,  the  first 
cavity  mode  can  be  calculated  from  Equation  2.  The  speed  of 
sound  in  the  water  column,  c0  is  first  computed  using  Equation 
5  presented  in  Chapter  III.  With  B  equal  to  2.18  x  109  N/nr  and 
Y3}  to  5.875  x  10iC  N/m2,  c0  is  found  to  be  846.06  m/s.  Q  is 
found  to  be  1.075  and  a  is  0.519.  With  these  values,  the 
cavity  resonance  frequency,  ojc,  is  7320.7  rad/s  or  fc  is  1165 
Hz.  The  relative  error  between  this  value  and  the  value 
obtained  from  the  modeling  (1250  Hz)  is  7%. 


D.  EFFECT  OF  ENCAPSULANT  MATERIAL  PROPERTIES  CHANGE 


The  ATILA  model,  without  the  mounting  plate,  was  used  to 
explore  the  effect  of  changing  the  material  properties  of  the 
encapsulant  materials ,  namely  the  polyurethane  around  the 
ceramic  and  the  neoprene  hoots. 

Figures  17  and  13  show  the  results  of  the  TVR  computed  for 
different  values  of  the  Young's  modulus  for  neoprene  and 
polyurethane,  respectively.  The  curves  are  plotted  over  the 
range  of  values  of  Young's  modulus  quoted  in  the 
manufacturer's  literature. 

TRANSMITTING  VOLTAGE  RESPONSE 

Varying  Polyurethane's  Young's  Modulus 


Neoprene:  Y=1E7  p^o.499  rho=14S2 


Figure  17:  Effects  of  Varying  Properties  of  Polyurethane 


TRANSMITTING  VOLTAGE  RESPONSE 

Varying  Neoprene's  Young's  Modulus 
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Figure  18:  Effect  of  Varying  Properties  of  Neoprene 


From  Figure  17,  it  can  be  seen  that  the  best  agreement 
with  the  measured  result  is  for  the  curve  with  Y  =  4.14  x  lO6 
N/m;  and  in  Figure  18,  it  was  the  curve  with  Y  =  1  x  10'  N/rrr. 
These  two  values  were  used  for  the  finite  element  in-water 
harmonic  analyses  presented  in  the  previous  sections. 

It  is  observed  from  the  two  figures  that  the  value  of 
Young's  modulus  for  polyurethane  has  much  greater  influence  on 
the  TVR  than  that  for  neoprene.  No  noticeable  change  in  ring 
resonance  frequency  was  observed  in  any  case.  The  transmitting 
voltage  response  was  found  to  increase  with  a  decrease  in  the 
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Young's  modulus  of  polyurethane.  It  was,  however,  the  reverse 
in  the  case  of  the  neoprene . 

E.  THE  ANALYTICAL  MODEL 
1.  The  Unloaded  Case 

The  expressions  developed  based  on  work  by  Hong- zhang 
Wang  [Ref.  2]  were  implemented  using  the  MATLAB  program.  The 
program  listing  is  provided  in  Appendix  3.  The  effective 
elastic  constant  values  obtained  for  the  unencapsulated, 
fiberglass  wrapped  ceramic  ring  assembly  of  the  finite  element 
model  were  used  for  the  analytical  model.  The  admittance 
circle  of  the  unencapsulated  fiberglass  wrapped  ring  is 
graphically  compared  to  the  measured  values  in  Figure  19. 


Figure  19:  Comparison  of  Analytical  with  Measured  Results  of 
Wrapped  Ring  In  Air 
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The  resonance  frequency  obtained  using  the  analytical 
method  -2420  Hz)  was  approximately  2  percent  lower  than  the 
measured  >2462  Hzs  and  the  value  obtained  using  the  finite 
element  code  ,2465  Hz)  . 

Wang's  approach  is  quicker  and  easier  to  use  compared 
-o  finite  element  modeling  as  only  the  circumferential 
compliance  sJ3  and  the  Poisson  ratio  are  required  for  the 
modeling.  It  is  especially  useful  when  only  measured 
admittance  circle  and  handbook  material  property  values  are 
available.  The  expressions  developed,  however,  do  not  take 
into  consideration  the  encapsulant  materials. 

Also,  Wang's  approach  could  not  be  applied  to  the 
piezoelectric  thin  ring  transducer  in  water,  as  it  was  not 
possible  to  obtain  expressions  for  the  radiation  impedances 
presented  to  the  inner  and  outer  surfaces  of  the  ring 
transducer. 
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VI.  CONCLUSIONS 


A  successful  two-dimensional  finite  element  model  of  a 
free- flooded,  tangentially  polarized  piezoelectric  ring  sonar 
transducer  was  built.  The  model  includes  internal  losses  of 
the  encapsulant  materials.  Although  the  model  includes  many 
necessary  simplifications  to  handle  the  problem  in  the 
available  MICROVAX  VMS  system,  they  were  designed  in  such  a 
way  that  the  transmitting  voltage  response  results  obtained 
from  the  in-water  harmonic  analyses  differs  by  only  2  dB  from 
the  measured  value  at  the  ring  resonance  frequency.  In  the 
region  of  the  cavity  mode,  the  model  including  the  mounting 
plates  has  better  agreement  except  for  a  dip  at  1500  Hz.  It 
was  able  to  accurately  predict  the  ring  resonance  frequency. 

The  model  built  (without  the  mounting  plate)  was  used  to 
investigate  the  effect  on  the  acoustical  properties  of  the 
transducer  due  to  a  change  in  the  properties  of  the 
encapsulant  materials  used  in  manufacturing.  It  was  observed 
that  the  Young's  modulus  of  the  polyurethane  has  a  significant 
influence  on  the  computed  transmitting  voltage  response. 
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For  the  analytical  model,  only  the  case  of  the  transducer 
in  air  was  developed.  The  resonance  frequency  was  found  to  be 
only  2  percent  lower  than  the  measured  value. 


VII.  SUGGESTIONS  FOR  FUTURE  WORK 

Suggestions  for  future  work  include: 

1-  Obtain  samples  of  the  actual  encapsulant  materials 
used  for  testing  to  determine  the  actual  material 
properties.  With  these  values,  the  model  can  then  be 
refined.  A  harmonic  analysis  should  then  be  performed 
and  the  transmitting  voltage  response  results  obtained 
compared  with  measured  values . 

2.  With  known  values  of  the  material  properties,  the 
effects  of  changes  in  encapsulant  materials  should 
again  be  investigated  using  the  in-water  model  with 
mounting  plate. 

3.  Finally,  it  would  be  very  helpful  and  interesting 
to  obtain  an  analytical  expression  for  the  radiation 
impedance  of  a  free- flooded  ring  transducer.  With  this 
and  the  results  of  Hong- zhang  Wang  [Ref.  2],  a  more 
complete  analytical  model  for  the  transducer  in  water 
can  be  produced.  The  results  obtained  with  this  model 
should  then  be  compared  to  those  obtained  using  ATILA 
and  the  measured  values , 

4.  Study  of  array  element  interaction  effects  should 
be  made  using  these  models  to  develop  techniques  for 
array  performance  modeling. 
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APPENDIX  A 


INPUT  DATA  FILE  FOR  HARMONIC  FINITE  ELEMENT  MODEL 


In-Water  Harmonic  Analysis  or  a  free-flooded  piezoelectric  ring 
sonar  transducer 

Prepared  By  MAJ  Tiong  Beng  Tay,  14  Feb  1993 


Navy  Type  I  ceramic 

Encapsulant  Material :  Polyurethane  &  Neoprene 
Tangentially  Polarized  but  modeled  as  radially-poled 
Radius  of  dipolar  damper  =  l  m 


RADIATION  DIPOLAR 
ANALYSIS  HARMONIC 
SKYLINE  COMPLEX 
PRECISION  DOUBLE 
CLASS  AXISYMME 
LCPDDC  =  7 

NLOAD  =40  *  The  number  of  loading  cases  must  be  less  than  or  equal 

*  to  100.  The  minimum  is  the  number  of  frequencies 

*  multiplied  by  2  because  we  are  analyzing  a  radiating 

*  transducer 

FREQUENCY 

0 . 225E+04  0 , 250E+04  0.275E+04  0.300E+04  0.325E+04  0.35E4  & 

0.375E4  0.4E4 

GEOMETRY  POLARIZA  CARTESIA 
1 

0 . 900E+Q2  0 . OQOE+OO  0.180E+03 

GEOMETRY 

2 

0 . 762E-03  0 . 000E+00 

3 

0 . 100E+01 


MATERIAL 

PZT4TA 


0. OOOE+OO 

0 . 000E+O0 

0 . 700E+04 

0. 000E+00 

0. 000E+00 

0. OOOE+OO 

& 

0 . 134E-10 

-0.443E-11 

-0 . 405E-11 

0. OOOE+OO 

0. 000E+00 

0. 000E+00 

& 

-0 .443E-11 

0 . 170E-10 

-0.443E-11 

0. 000E+00 

0. 000E+00 

0. OOOE+OO 

& 

-0  .4QSE- 11 

-0.443E-11 

0.134E-10 

0. OOOE+OO 

0. OOOE+OO 

0. 000E+00 

& 

0 .000E+00 

0. 000E+00 

0. 000E+00 

0.390E-10 

0. 000E+00 

0. 000E+00 

& 

0 . OQOE+OO 

0. OOOE+OO 

0. 000E+00 

0. 000E+00 

0 . 357E-10 

0. 000E+00 

& 

0 . OOOE+QO 

0. 000E+0O 

0. 000E+00 

0. OOOE+OO 

0. 000E+00 

0.390E-10 

Sc 

0 . 000E+00 

0 .OOOE+OO 

0. OOOE+OO 

0 .000E+00 

0.525E-09 

0. OOOE+OO 

Sc 

0 . 000E+00 

0. OOOE+OO 

0. 000E+00 

0 .525E-09 

0. 000E+00 

0. OOOE+OO 

Sc 

-0 . 135E-09 

0 . 300E-09 

-0.135E-09 

0 .000E+00 

0. 000E+00 

0. 000E+00 

Sc 

0.624E-08 

0. 000E+00 

0. 000E+00 

0. OOOE+OO 

0. 000E+00 

0 .000E+00 

Sc 

0 . OOOE+OO 

0.624E-08 

0. OOOE+OO 

0. 000E+00 

0. OOOE+OO 

0. 000E+00 

Sc 

0 . Q00E+0Q 

0. 000E+00 

0.570E-08 

0.  OOOE+OO 

0. 000E+00 

0. OOOE+OO 

FGLR 

0 . 587E+11 

0.280E+00 

0 . 254E+04 

P1S90 

0.100E+09 

0 . 490E+00 

0.110E+04 

0.  OOOE+OO 

0. 000E+00 

0.400B+08 
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NEOPRENE 


0  . 

100E+08  0 . 499E+00 

0 . 148E+04 

HRS 

n 

195E+1 

2  0 . 333E+00 

0 . 770E+04 

ENTRAN 

0  . 

X38E+08  0 . 330E+00 

0 . 170E+04 

NODES 

* 

1 

k 

-0 . 86360E-01 

0 . 18339E+00 

nr 

4 

* 

-0 .86360E-01 

0 . 18656E+00 

3 

★ 

-0 .86360E-01 

0 . 18974E+00 

k 

4 

★ 

-  0 . 75565E- 01 

0 . 18339E+00 

k 

5 

k 

-  0 . 75565E-01 

0 . 18974E+00 

★ 

6 

★ 

-  0 . 64770E- 01 
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k 

444 

*  0 . 86360E- 01 

0 , 11500E+00 

0 . OOOOOE+OO 

k 

445 

*  -  0 .  86360E-01 

0 . 11500E+00 

0 . OOOOOE+OO 

k 

446 

*  -0.10795E+00 

0.98290E-01 

0. OOOOOE+OO 

k 

447 

*  0 . 10795E+00 

0 , 98290E-01 

0. OOOOOE+OO 

k 

448 

*  0 . 14605E+00 

0. OOOOOE+OO 

0. OOOOOE+OO 

k 

449 

*  - 0 . 14605E+00 

0.  OOOOOE+OO 

0  .OOOOOE+OO 

k 

450 

*  0 . 64770E- 01 

0 . 13170E+00 

0. OOOOOE+OO 

k 

451 

*  -  0 . 64770E-01 

0  - 13170E+00 

0. OOOOOE+OO 

k 

452 

*  0 . 14605E+00 

0.24576E-01 

0. OOOOOE+OO 

k 

453 

*  -0.14605E+00 

0 , 24576E-01 

0. OOOOOE+OO 

k 

454 

*  0 . OOOOOE+OO 

0 . 14840E+00 

0. OOOOOE+OO 

k 

455 

*  0.95250E-01 

0 . 11500E+00 

0 .OOOOOE+OO 

k 

456 

*  -0.95250E-01 

0 .11500E+00 

0. OOOOOE+OO 

k 

457 

*  0.21590E-01 

0 . 14840E+00 

0 . 000008+00 

k 

458 

*  -Q.21590E-01 

0 . 14840E+00 

0.  OOOOOE+OO 

k 

459 

*  -0.11430E+00 

0.98290E-01 

0. OOOOOE+OO 

k 

460 

*  0  - 11430E+00 

0.98290E-01 

0. OOOOOE+OO 

k 

461 

*  0 . 75565E-01 

0 . 13170E+Q0 

0. OOOOOE+OO 

k 

462 

*  -0 . 75565E-01 

0.13170E+00 

0. OOOOOE+OO 

k 

463 

*  * 0 . 15240E+00 

0. OOOOOE+OO 

0 .OOOOOE+OO 

k 

464 

*  0 . 14605E+00 

0.49150E-01 

0. OOOOOE+OO 

k 

465 

*  -0.14605E+00 

0.49150E-01 

0 .OOOOOE+OO 

k 

466 

*  0.43180E-01 

0.14840E+00 

0. OOOOOE+OO 

k 

467 

*  -0.43180E-01 

0 . 14840E+00 

0. OOOOOE+OO 

k 

468 

*  -0 . 12065E+00 

0 . 98290E-01 

0. OOOOOE+OO 

k 

469 

*  0.12065E+00 

0 , 98290E- 01 

0. OOOOOE+OO 

k 

470 

*  0 . 86360E-01 

0 . 13170E+00 

0. OOOOOE+OO 

k 

471 

*  -0 . 86360E-01 

0 . 13170E+00 

0. OOOOOE+OO 

k 

472 

*  0 . 10795E+00 

0 . 115QQE+00 

0. OOOOOE+OO 

k 

473 

*  '0 . 10795E+00 

0 . 11500E+00 

0. OOOOOE+OO 

k 

474 

*  '0 . 15875E+00 

0. OOOOOE+OO 

0. OOOOOE+OO 

k 

475 

*  0.90805E-01 

0 . 13170E+00 

0. OOOOOE+OO 

* 

476 

*  - 0 . 90805E-01 

0 . 13170E+00 

0. OOOOOE+OO 

k 

477 

*  -0 . 15240E+00 

0 . 49150E-01 

0 . OOOOOB+OO 

k 

478 

*  -0 . 15875E+00 

0 . 24576E-01 

0. OOOOOE+OO 

k 

479 

*  0 . 64770E- 01 

0 . 14840E+00 

0. OOOOOB+OO 

k 

480 

*  -0.64770B-01 

0 . 14840E+00 

0. OOOOOE+OO 

* 

481 

*  0.95250H-01 

0 . 13170E+00 

0. OOOOOB+OO 

k 

482 

*  -0.95250B-01 

0 . 13170B+00 

0. OOOOOB+OO 

Jr 

483 

*  0 . 14605E+00 

0.73721H-01 

0. OOOOOB+OO 

k 

484 

*  -0.14605E+00 

0.73721E-01 

0. OOOOOE+OO 

k 

485 

*  0. 00000E+00 

0 . 16510E+00 

0.  OOOOOB+OO 

k 

486 

*  Q.10795B-Q1 

0.16510B+00 

0.  OOOOOE+OO 

k 

487 

*  -0.10795E-01 

0.16510B+00 

0. OOOOOB+OO 
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*  488  *  - 0 . 13350E+GQ  0.98290E-01  O.OOOOOE+OO 

'  489  *  0 . 13350E+00  0.98290E-01  O.OOOOOE+OO 

*  490  *  -  0 . 15875E+00  0.49150E-01  O.OOOOOE+OO 

*  491  *  0 . 10160E+00  0 . 13170E+00  O.OOOOOE+OO 

*  492  *  -0 . 1Q160E+00  0.13170E+00  O.OOOOOE+OO 

*  493  *  0 . 21590E- 01  0.16510E+00  O.OOOOOE+OO 

*  494  *  -  0 . 21590E  -  01  0.16510E+00  O.OOOOOE+OO 

*  495  *  0 . 12065E+00  0.11500E+00  O.OOOOOE+OO 

*  496  *  -0 . 12065E+00  0 . 11500E+00  O.OOOOOE+OO 

*  497  *  0.32385E-01  0.16510E+00  O.OOOOOE+OO 

*  498  *  -0 . 32385E-01  0.16510E+00  O.OOOOOE+OO 

*  499  *  0 . 10795E+00  0.13170E+00  O.OOOOOE+OO 

*  500  *  -  0 . 10795E+00  0.13170E+00  O.OOOOOE+OO 

*  501  *  0.43180E-01  0.16510E+00  O.OOOOOE+OO 

*  502  *  -0.43180E-01  0.16510E+00  O.OOOOOE+OO 

*  503  *  0 . 86360E -  01  0.14840E+00  O.OOOOOE+OO 

*  504  *  -  0 . 86360E- 01  0.14840E+00  O.OOOOOE+OO 

*  505  *  0 . 53975E- 01  0.16510E+00  O.OOOOOE+OO 

*  506  *  -  0 . 53975E- 01  0.16510E+00  O.OOOOOE+OO 

*  507  *  0 . 11430E+00  0.13170E+00  O.OOOOOE+OO 

*  508  *  -  0 . 11430E+00  0.13170E+00  O.OOOOOE+OO 

*  509  *  -  0 . 15875E+00  0.73721E-01  O.OOOOOE+OO 

*  510  *  0 . 14605E+00  0.98290E-01  O.OOOOOE  +  OO 

*  511  *  - 0 . 14605E+00  0.98290E-01  O.OOOOOE+OO 

*  512  *  0.95250E-01  0.14840E+00  O.OOOOOE+OO 

*  513  *  -  0 . 95250E- 01  0.14840E+00  O.OOOOOE+OO 

*  514  *  0.64770E-01  0.16510E+00  O.OOOOOE+OO 

*  515  *  - 0 . 64  770E- 01  0.16510E+00  O.OOOOOE+OO 

*  516  *  0 . 12065E+00  0.13170E+00  O.OOOOOE+OO 

*  517  *  -0 . 12065E+00  0.13170E+00  O.OOOOOE+OO 

*  518  *  -  0 . 15240E+00  0.98290E-01  O.OOOOOE+OO 

*  519  *  0 . 75565E- 01  0.16510E+00  O.OOOOOE+OO 

*  520  *  -  0 , 75565E-01  0.16510E+00  O.OOOOOE+OO 

*  521  *  0 . 10795E+00  0.14840E+00  O.OOOOOE+OO 

*  522  *  - 0 . 10795E+00  0.14840E+00  O.OOOOOE+OO 

*  523  *  0 . 14605E+00  0.11500E+00  O.OOOOOE+OO 

*  524  *  - 0 . 14605E+00  0.11500E+00  O.OOOOOE+OO 

*  525  *  0 . 86360E-01  0.16510E+00  O.OOOOOE+OO 

*  526  *  -0 . 86360E-01  0.16510E+00  O.OOOOOE+OO 

*  527  *  -0 . 15875E+00  0.98290E-01  O.OOOOOE+OO 

*  528  *  0 . 13350E+00  0.13170E+00  O.OOOOOE+OO 

*  529  *  - 0 . 13350E+00  0.13170E+00  O.OOOOOE+OO 

*  530  *  0.90805E-01  0.16510E+00  O.OOOOOE+OO 

*  531  *  -0.90805E-01  0.16510E+00  O.OOOOOE+OO 

*  532  *  0 . 18875E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  533  *  0.95250E-01  0.16510E+00  O.OOOOOE+OO 

*  534  *  -  0 . 95250E- 01  0.16510E+00  O.OOOOOE+OO 

*  535  *  0 . 12065E+00  0.14840E+00  O.OOOOOE+OO 

*  536  *  -0 . 12065E+00  0.14840E+00  O.OOOOOE+OO 

*  537  *  0 . 10160E+00  0.16510E+00  O.OOOOOE+OO 

*  538  *  -0 . 1Q16QE+00  0.16510E+00  O.OOOOOE+OO 

*  539  *  0 . 18869E+00  0.4915QE-01  O.OOOOOE+OO 

*  540  *  -0 . 15875E+00  0.11500E+00  O.OOOOOE+OO 

*  541  *  0 . 14605E+00  0.13170E+00  O.OOOOOE+OO 

*  542  *  -0 . 14605E+00  0.13170E+00  O.OOOOOE+OO 

*  543  *  0 . 10795E+00  0.16510E+00  O.OOOOOE+OO 

*  544  *  - 0 . 10795E+00  0.16510E+00  O.OOOOOE+OO 

*  545  *  0 . 11430E+00  0.16510E+00  O.OOOOOE+OO 

*  546  *  -0 . 11430E+00  0.16510E+00  O.OOOOOE+OO 

*  547  +  -0.20081E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  548  *  0.12065E+00  0.I6510E+00  O.OOOOOE+OO 

*  549  *  -0 . 12065E+00  0.16510E+00  O.OOOOOE+OO 
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550 

★ 

0  .OOOOOE+OO 

* 

551 

k 

-  0 . 10795E- 01 

* 

552 

k 

0 . 10795E- 01 

★ 

553 

k 

-0 . 21590E-01 

★ 

554 

k 

0 . 21590E- 01 

* 

555 

k 

-0 . 15875E+00 

556 

k 

-0 .20075E+00 

★ 

557 

k 

-  0 . 32385E -  01 

★ 

558 

k 

0.32385E-01 

* 

559 

k 

0 . 14605E+00 

* 

560 

k 

-0.14605E+00 

★ 

561 

k 

-0 .43180E-01 

★ 

562 

k 

0.43180E-01 

★ 

563 

k 

-  0 . 53975E- 01 

★ 

564 

k 

0.53975E-01 

* 

565 

k 

0 . 13350E+00 

★ 

566 

k 

-0 .13350E+00 

★ 

567 

k 

0 . 18851E+00 

it 

568 

k 

-0 . 64770E-01 

k 

569 

k 

0 . 64770E- 01 

k 

570 

k 

-  0 . 15875E+00 

it 

571 

k 

-0 . 75565E-01 

it 

572 

k 

0 . 75565E-01 

k 

573 

k 

0 . 14605E+00 

k 

574 

k 

-0 . 14605E+00 

k 

575 

k 

-0.86360E-01 

* 

576 

k 

0 . 863fa0E-0x 

★ 

577 

k 

-0.20057E+00 

* 

578 

k 

0 . 14605E+00 

k 

579 

k 

-0.90805E-01 

k 

580 

k 

0.90805E-01 

k 

581 

k 

-0.95250E-01 

k 

582 

k 

0.95250E-01 

k 

583 

k 

0 . 14605E+00 

k 

584 

k 

-0.10160E+00 

k 

585 

k 

0 . 10160E+00 

k 

586 

k 

-0 . 15875E+00 

k 

587 

k 

0 . 18831E+00 

k 

588 

k 

0 . 14605E+00 

k 

589 

k 

-0.10795E+00 

k 

590 

k 

0 . 10795E+00 

k 

591 

k 

0.23144E+00 

k 

592 

k 

-0.15875E+00 

k 

593 

k 

0.23141E+00 

k 

594 

k 

0.14605E+00 

k 

595 

k 

-0 . 11430E+00 

k 

596 

k 

0 . 11430E+00 

k 

597 

k 

-  0 . 15875E+00 

k 

598 

k 

0 .23132B+00 

k 

599 

k 

0 . 14605E+00 

k 

600 

k 

0.12065B+00 

k 

601 

k 

-0 , 12065E+00 

k 

602 

k 

-0 . 15875E+00 

k 

603 

k 

0 . 14605E+00 

k 

604 

k 

0 . 14605E+00 

k 

605 

k 

-0.15875E+00 

k 

606 

k 

0 . 23117E+00 

k 

607 

k 

-0 . 24288B+00 

k 

608 

k 

0 . 14605B+00 

k 

609 

k 

-0 . 24284E+00 

k 

610 

k 

-0 . 13350B+00 

k 

611 

k 

0.13350E+00 

0 . 20472E+0Q 
0 . 20472E+00 
0 . 20472E+00 
0.20472E+00 
0 .20472E+00 
0 . 13170E+00 
0 . 49150E-01 
0 . 20472E+00 
0 . 20472E+00 
0 . 14840E+00 
0.14840E+00 
0 . 20472E+00 
0 . 20472E+00 
0 . 20472E+00 
0 .20472E+00 
0 .16510E+00 
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0 . 14840E+00 
0.20472E+00 
0 . 20472E+00 
0.16510E+00 
0.16510E+00 
0.20472E+00 
0 . 20472E+00 
0.98290E-01 
0 . 16942E+00 
0 . 20472E+00 
0 . 2Q472E+00 
0 . 20472E+00 
0 . 20472E+00 
0.17374B+00 
0.20472E+00 
0 . 20472E+00 
0 . 16510E+00 
0 . 13170E+00 
0 .17856E+00 
0 . 20472E+00 
0 .20472E+00 
0 . OOOOOE+OO 
0 . 16942E+00 
0.24576E-01 
0 . 18339E+00 
0 . 20472E+00 
0 . 20472E+00 
0 . 17374E+00 
0.49150E-01 
0 . 18656E+00 
0 . 20472E+QQ 
0 . 20472E+00 
0 . 17856E+00 
0 . 18974E+00 
0 . X9202E+00 
0 . 18339E+00 
0 . 73721E-01 
0 . OOOOOE+OO 
0 . X9430E+00 
0.24576B-01 
0 . 20472E+00 
0.20472B+00 
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0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
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612 

ie 

0  .OOOOOE+OO 

0 .24449E+00 

613 

-  0 . 15875E+00 

0 . 18656E+00 

614 

* 

0 . 22142E -  01 

0 .24446E+00 

615 

★ 

-0 .22142E-01 

0 .24446E+00 

616 

w 

0 . 14605E  +  00 

0 . 19951E+00 

617 

-0 . 15875E+00 

0 . 18974E+00 

618 

* 

-0 . 24275E+00 

0 .49150E-01 

619 

■*r 

0 .44285E-01 

0 .24438E+00 

620 

•* 

-0 .44285E-01 

0 .24438E+00 

621 

* 

-0 . 15875E+00 

0 . 19202E+00 

622 

ft 

0 . 18806E+00 

0 . 16510E+00 

623 

*■ 

-0 .15875E+00 

0 .19430E+00 

624 

nr 

0 .23096E+00 

0 .98290E-01 

625 

★ 

0 .14605E+00 

0 . 20472E+00 

62  6 

w 

-0 . 14605E+00 

0 . 20472E+00 

627 

★ 

0 . 66427E- 01 

0 .24425E+00 

628 

★ 

-0 .66427E-01 

0 .24425E+00 

629 

★ 

-  0 . 24260E+00 

0 . 73721E-01 

630 

★ 

-0 .15875E+00 

0 . 19951E+00 

631 

-  0 . 15240E+00 

0 .20472E+00 

632 

* 

0 .23078E+00 

0 . 11500E+00 

633 

★ 

0 . 18704E+00 

0 . 17799E+00 

634 

* 

-0 .1587SE+00 

Q . 2Q472E+00 

635 

* 

-0 .20013E+00 

0 . 16510E+00 

636 

ft 

0 . 88569E-01 

0 .24406E+00 

637 

★ 

-  0 . 88569E- 01 

0 .24406E+00 

638 

* 

-  0 . 24239E+00 

0 .98290E-01 

639 

* 

-0 .10088E+00 

0 . 24340E+00 

640 

★ 

0 . 10301E+00 

0 .24290E+00 

641 

★ 

0.23057E+00 

0 . 13170E+00 

642 

★ 

-  0 . 19911E+00 

0 . 17799E+00 

643 

★ 

0 . 18406E+00 

0 . 19535E+00 

644 

★ 

-  0 . 11549E+00 

0 . 24278E+00 

645 

★ 

0 . 12088E+00 

0 . 24101E+00 

646 

•k 

0 . 23034E+00 

0 . 14840E+00 

647 

★ 

0 . 27414E+00 

0. OOOOOE+OO 

648 

★ 
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0 . 24101E+00 

649 

* 

-0 .24200E+00 

0 . 13170E+00 
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★ 
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* 
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★ 
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★ 
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655 

★ 

0 . 23007E+00 

0 . 16510E+Q0 
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★ 

-0 . 19177E+00 

0 . 20893E+00 
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* 
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★ 
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★ 
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★ 
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* 
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★ 
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0 . 28420E+00 
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* 
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★ 

-0 . 34042E-01 
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★ 
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* 
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0.28404E+00 
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* 

-  0 . 45389E- 01 
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★ 
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* 
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0. OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0 .OOOOOE+OO 
0 .OOOOOE+OO 


674  *  0.27342E+00  G.9829QE-Q1  O.OOOOOE+OO 

575  *  0 . 68084E-01  0.28377E+00  O.OOOOOE+OO 

676  *  -0 . 68084E-01  0.28377E+00  O.OOOOOE+OO 

677  *  0 . 22804E+00  0 . 18225E+00  O.OOOOOE+OO 

678  *  -  0 . 24151E+00  0 . 16510E+0Q  O.OOOOOE+OO 

679  *  Q.79431E-0X  0.28360E+00  O.OOOOOE+OO 

680  *  -Q.79431E-01  0.28360E+00  O.OOOOOE+OO 

681  *  - 0 . 179 85E+00  0.23519E+00  O.OOOOOE+OO 

682  *  -  0 . 24060E+00  0.17370E+00  O.OOOOOE+OO 

683  *  0 . 90779E -  01  0.28340E+00  O.OOOOOE+OO 

684  *  -0.90779E-01  0.28340E+00  O.OOOOOE+OO 

685  *  0 . 22542E+00  0.19491E+00  O.OOOOOE+OO 

686  *  -  0 . 98660E -  01  0.28281E+00  O.OOOOOE+OO 

687  *  0 . 10081E+Q0  0.28241E+00  O.OOOOOE+OO 

688  *  - 0 . 28421E+00  0.98290E-01  O.OOOOOE+OO 

689  *  -  0 . 2394  7E+00  0.18225E+00  O.OOOOOE+OO 

690  *  - 0 . 10651E+00  0.28207E+00  O.OOOOOE+OO 

691  *  0 . 11076E+00  0.28107E+00  O.OOOOOE+OO 

692  *  0 . 27284E+00  0.13171E+00  O.OOOOOE+OO 

693  *  0 . 22208E+00  0.20731E+00  O.OOOOOE+OO 

694  *  - 0 . 11478E+00  0.28149E+00  O.OOOOOE+OO 

695  *  0 . 12234E+00  0.27936E+00  O.OOOOOE+OO 

696  *  - 0 . 12304E+00  0.28084E+00  Q.OOOOOE+OO 

697  +  -  0 . 23685E+00  0.19491E+00  O.OOOOOE+OO 

698  *  0 . 13380E+00  0.27729E+00  O.OOOOOE+OO 

699  *  0.21805E+OO  0.21791E+00  O.OOOOOE+OO 

700  *  -0 . 13419E+Q0  0.27922E+00  O.OOOOOE+OO 

701  *  0.14513E+00  0 . 27488E+00  O.OOOOOE+OO 

702  *  -  0 . 23350E+00  0.20731E+00  O.OOOOOE+OO 

703  *  0 . 21336E+00  0.22812E+00  O.OOOOOE+OO 

704  *  -0 . 14524E+00  0.27729E+00  O.OOOOOE+OO 

705  *  0 . 15630E+00  0.27213E+00  O.COOOOE+OO 

706  *  0 . 20805E+00  0.23710E+00  O.OOOOOE+OO 

707  +  -  0 . 22948E+00  0.21791E+00  O.OOOOOE+OO 

708  *  0 . 31684E+00  O.OOOOOE+OO  O.OOOOOE+OO 

709  *  0.31678E+00  0.24576E-01  O.OOOOOE+OO 

710  *  0 . 20216E+00  0.24558E+00  O.OOOOOE+OO 

711  *  0 . 27209E+00  0.16510E+00  O.OOOOOE+OO 

712  *  -0 . 16241E+00  0 . 27488E+00  O.OOOOOE+OO 

713  *  0 . 17314E+00  0.26906E+00  O.OOOOOE+OO 

714  *  -  0 . 22479E+00  0.22812E+00  O.OOOOOE+OO 

715  *  0 . 31660E+00  0.49150E-01  O.OOOOOE+OO 

716  *  0 . 19606E+00  0.25588E+00  O.OOOOOE+OO 

717  *  - 0 . 21948E+00  0.23710E+00  O.OOOOOE+OO 

718  *  O.OOOOOE+OO  0.32402E+00  O.OOOOOE+OO 

719  *  0.23247E-01  0.32394E+00  O.OOOOOE+OO 

720  *  -Q.23247E-01  0.32394E+Q0  O.OOOOOE+OO 

721  *  0 . 31630E+00  0.73721E-01  O.OOOOOE+OO 

722  *  -0  -21359E+00  0.24558E+00  O.OOOOOE+OO 

723  *  -0.17916E+00  0.27213E+00  O.OOOOOE+OO 

724  *  0 . 18951E+00  0.26566E+00  O.OOOOOE+OO 

725  *  -0 . 32700E+00  O.OOOOOE+OO  O.OOOOOE+OO 

726  *  0.46494E-01  0.32369E+00  O.OOOOOE+OO 

727  *  -0 . 46494E-01  0.32369E+00  O.OOOOOE+OO 

728  *  0.26903E+00  0.18650E+00  O.OOOOOB+OO 

729  *  -0 . 28288E+00  0.16510B+00  O.OOOOOB+OO 

730  *  -Q.32694E+00  0.24576E-01  O.OOOOOE+OO 

731  *  - 0 . 20749E+00  0.255888+00  O.OOOOOB+OO 

732  *  -0.19015E+00  0.26906E+00  O.OOOOOE+OO 

733  *  -0.32676E+00  0.49149E-01  O.OOOOOE+OO 

734  *  0 . 69741E-01  0.32329B+00  O.OOOOOB+OO 

735  *  -0 . S9741E-01  0.32329E+00  O.OOOOOE+OO 


-  ”36  *  0 . 31587E+00  0.98290E-01  O.OOOOOE+OO 

*  737  *  -  0 . 20094E+00  0 . 26566E+00  O.OOOOOE+OO 

*  738  *  -  0 . 32646E+00  0.73720E-01  j.GOOOOE+OO 

'  739  *  0 . 3 1551E+00  0 . X1500E+00  O.OOOOOE  +  OO 

*  ”40  *  0.929S8E-01  0.32273E+00  O.OOOOOE+OO 

*  ”41  *  -  0 . 92988E-01  0.32273E+00  O.OOOOOE+OO 

'  ”42  ~  - 0 . 27983E+00  0.18650E+00  C.OOOOOE+OO 

*  743  *  - 0 . 11215E+00  0 . 32074E+00  O.OOOOOE+OO 

*  744  *  0 . 26009E+00  0.21927E+00  O.OOOOOE+OO 

-  745  *  -  0 . 32603E+00  0.98290E-01  O.OOOOOE  +  OO 

*  746  *  0.11852E+00  0.31924E+00  O.OOOOOE+OO 

'  747  *  0 . 3 1510E+00  0.13171E+00  O.OOOOOE+OO 

*  748  *  -  0 . 13058E+00  0.31890E+00  O.OOOOOE+OO 

*  749  *  0 . 14673E+00  0.31358E+00  O.OOOOOE+OO 

*  750  *  0 . 31463E+00  0.14841E+00  O.OOOOOE+OO 

*  751  *  -0 . 27088E+00  0.21927E+00  O.OOOOOE+OO 

*  752  *  0 . 24701E+00  0.24731E+00  O.OOOOOE+OO 

*  753  *  - 0 . 32526E+00  0.13171E+00  O.OOOOOE+OO 

*  754  *  -  0 . 15753E+00  0.31358E+00  O.OOOOOE+OO 

*  755  *  0 . 17413E+00  0.30584E+00  O.OOOOOE+OO 

*  756  *  0 . 31410E+00  0.16510E+00  O.OOOOOE+OO 

*  757  *  0 . 23021E+00  0.27122E+00  O.OOOOOE+OO 

'  75g  *  -  0 . 2S781E+00  0. 24731E+00  O.OOOOOE+OO 

»  ”59  *  0 . 31229E+00  0.17798E+00  O.OOOOOE+OO 

*  760  *  0 . 35954E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  761  *  0 . 35923E+00  0.49150E-01  O.OOOOOE+OO 

*  762  *  - 0 . 24100E+00  0.27122E+00  O.OOOOOE+OO 

*  763  *  -0 . 19572E+00  0.30584E+00  O.OOOOOE+OO 

*  764  *  0 . 21125E+00  0.29613E+00  O.OOOOOE+OO 

*  765  *  O.OOOOOE+OO  0.36378E+00  O.OOOOOE+OO 

*  766  *  - 0 . 32426E+00  0.16510E+00  O.OOOOOE+OO 

*  767  *  0 . 11900E- 01  0.36375E+00  O.OOOOOE+OO 

*  768  *  -0 . 11900E-01  0.36375E+00  O.OOOOOE+OO 

*  769  *  0 . 31003E+00  0.19076E+00  O.OOOOOE+OO 

*  770  *  0.23799E-01  0.36367E+00  O.OOOOOE+OO 

*  771  *  -0.23799E-01  0.36367E+00  O.OOOOOE+OO 

*  772  *  0 . 35699E-01  0.36354E+00  O.OOOOOE+OO 

*  773  *  -0.35699E-01  0.36354E+00  O.OOOOOE+OO 

*  774  *  0.47599E-01  0.36335E+00  O.OOOOOE+OO 

*  775  *  - 0 . 47599E - 01  0.36335E+00  O.OOOOOE+OO 

*  776  *  0 . 59498E- 01  0.36311E+00  O.OOOOOE+OO 

*  777  *  - 0 . 59498E- 01  0.36311E+00  O.OOOOOE+OO 

*  778  *  - 0 . 32245E+00  0.17798E+00  O.OOOOOE+OO 

*  779  *  -0.36906E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  780  *  0.71398E-01  0.36282E+00  O.OOOOOE+OO 

*  781  *  -0 . 71398E-01  0.36282E+00  O.OOOOOE+OO 

*  782  *  - 0 . 222Q4E+00  0.29613E+00  O.OOOOOE+OO 

*  783  *  0 . 30478E+00  0.21125E+00  O.OOOOOE+OO 

*  784  *  0 . 35833E+00  0.98290E-01  O.OOOOOE+OO 

*  785  *  0 . 83298E-01  0.36247E+00  O.OOOOOE+OO 

*  786  *  -0 . 83298E-01  0.36247E+00  O.OOOOOE+OO 

*  787  *  -0.36876E+00  0.49149E-01  O.OOOOOE+OO 

*  788  *  - 0 . 32019E+00  0.19076E+00  O.OOOOOE+OO 

*  789  *  0 . 95197E-01  0.36207E+00  O.OOOOOE+OO 

*  790  *  -0.95197E-01  0.36207E+00  O.OOOOOE+OO 

*  791  *  -0 . 10652E+00  0.36089E+00  O.OOOOOE+OO 

*  792  *  0 . 11081E+00  0.36011E+00  O.OOOOOE+OO 

*  793  *  0 . 29810E+00  0.23123E+00  O.OOOOOE+OO 

*  794  *  -0 . 11778E+00  0.35941E+00  O.OOOOOE+OO 

*  795  *  0.12628E+00  0.35741E+00  O.OOOOOE+OO 

*  796  *  -0 . 31494E+00  0.21126E+00  O.OOOOOE+OO 

*  797  *  - 0 . 12797E+00  0.35826E+00  O.OOOOOE+OO 
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HP 

^93 

*  -  0 . 36785E+00 

0 .98290E-01 

0  .  OOOOOE+OO 

■>r 

799 

*  0 . 35736E+00 

0 . 13171E+00 

0  .  OOOOOE+OO 

k 

aoo 

*  0 . 14308E+00 

0 . 35399E+00 

0 . OOOOOE+OO 

k 

801 

*  0 . 29004E+00 

0 . 24925E+00 

0 .OOOOOE+OO 

* 

302 

*  -0 . 13812E+00 

0.35697E+00 

0  .  OOOOOE+OO 

* 

803 

*  0 . 15966E+00 

0 . 34986E+00 

0. OOOOOE+OO 

★ 

804 

*  -  0 . 30826E+00 

0 . 23123E+00 

0 .  OOOOOE+OO 

★ 

805 

*  -  0 . 15407E+00 

0 . 35372E+00 

0 . OOOOOE+OO 

* 

806 

*  0 . 28067E+00 

0 . 26651E+00 

0  .OOOOOE+OO 

★ 

307 

*  0 . 17596E+00 

0 . 34504E+00 

0 .OOOOOE+OO 

★ 

808 

*  - 0 . 16982E+00 

0 . 34986E+00 

0 .OOOOOE+OO 

★ 

809 

*  0 . 19195E+00 

0 . 33954E+00 

0  .  OOOOOE+OO 

★ 

810 

*  -  0 . 30020E+00 

0 .24925E+00 

0  .OOOOOE+OO 

* 

311 

*  0 . 27004E+00 

0 .28218E+00 

0  .  OOOOOE+OO 

★ 

812 

*  0 . 35611E+00 

0 .16510E+00 

0  .  OOOOOE+OO 

k 

813 

*  0 . 25826E+00 

0 .29686E+00 

0. OOOOOE+OO 

★ 

814 

*  - 0 . 29Q83E+00 

0 .26651E+00 

0. OOOOOE+OO 

k 

815 

*  -  0 . 19132E+00 

0 .34504E+00 

0. OOOOOE+OO 

★ 

816 

*  0 .  21278E+00 

0 . 33339E+00 

0  .OOOOOE+OO 

* 

817 

*  0 . 24607E+00 

0 . 31226E+00 

0 .  OOOOOE+OO 

■k 

818 

*  -  0 . 28020E+00 

0 . 28218E+00 

0 .  OOOOOE+OO 

-k 

819 

*  - 0 . 26842E+00 

0 . 29686E+00 

0. OOOOOE+OO 

k 

820 

*  - 0 . 21227E+00 

0 . 33954E+00 

0. OOOOOE+OO 

k 

821 

*  -  0 . 36564E+00 

0 . 16510E+00 

0. OOOOOE+OO 

k 

822 

*  0 . 23298E+00 

0 . 32660E+00 

0.  OOOOOE+OO 

k 

823 

*  0 . 35102E+00 

0 . 19501E+00 

0. OOOOOE+OO 

★ 

824 

*  0 . 4  0224E+00 

0 . OOOOOE+QO 

0. OOOOOE+OO 

★ 

825 

*  0 . 40214E+00 

0 .24576E-01 

0. OOOOOE+OO 

k 

826 

*  0  .  QOOOQE+OO 

0.40354E+00 

0.  OOOOOE+OO 

k 

827 

*  -0 . 22790E+Q0 

0 . 33339E+00 

0. OOOOOE+OO 

k 

828 

*  -0.25622E+00 

0 . 31226E+00 

0. OOOOOE+OO 

k 

829 

*  0 . 24352E-01 

0.40341E+00 

0. OOOOOE+OO 

k 

830 

*  -0.24352E-01 

0.40341E+00 

0. OOOOOE+OO 

k 

831 

*  0 . 40187E+00 

0 . 49150E-01 

0. OOOOOE+OO 

k 

832 

*  0.48703E-01 

0 . 40301E+00 

0. OOOOOE+OO 

k 

833 

*  -0.48703E-01 

0 . 40301E+00 

0 . OOOOOE+OO 

k 

834 

*  - 0 . 24314E+00 

0.32660E+00 

0. OOOOOE+OO 

k 

835 

*  0 . 40142E+00 

0 . 73721E-01 

0 .OOOOOE+OO 

k 

836 

*  0 . 73055E-01 

0 . 40234E+00 

0 . OOOOOE+OO 

k 

837 

*  -  0 . 73055E- 01 

0 . 40234E+00 

0 . OOOOOE+OO 

k 

838 

*  - 0 .  36055E+00 

0 . 19501E+00 

0. OOOOOE+OO 

k 

839 

*  -  0 .  41113E+00 

0 .OOOOOE+OO 

0 .OOOOOE+OO 

k 

840 

*  -0.41103E+00 

0.24575E-01 

0 -OOOOOE+OO 

k 

841 

*  0.40078E+00 

Q.98290E-01 

0. OOOOOE+OO 

k 

842 

*  0 . 97406E- 01 

0 .40140E+00 

0 -OOOOOE+OO 

k 

843 

*  -0.97406E-01 

0 . 40140E+00 

0. OOOOOE+OO 

k 

844 

*  -0.41076E+00 

0 . 49149E-01 

0 . OOOOOE+OO 

k 

845 

*  0 . 33612E+00 

0 . 24319E+00 

0. OOOOOE+OO 

k 

846 

*  0 . 40024E+00 

0 . 11501E+00 

0. OOOOOE+OO 

k 

847 

*  -0 . 12341E+00 

0.39808E+00 

0. OOOOOE+OO 

k 

848 

*  -0.41031E+00 

0.73720E-01 

0. OOOOOE+OO 

k 

849 

*  0 . 13403E+00 

0 . 39558E+00 

0. OOOOOE+OO 

k 

850 

*  0 . 39962E+00 

0 . 13171E+00 

0- OOOOOE+OO 

k 

851 

*  -Q.14567E+00 

0 . 39503E+00 

0. OOOOOE+OO 

k 

852 

*  -0.40967E+00 

0 . 98290E-01 

0 .OOOOOB+OO 

k 

853 

*  -0 . 34564E+00 

0 . 24319E+00 

0. OOOOOE+OO 

k 

854 

*  0 , 17259E+00 

0 . 38615E+00 

0. OOOOOB+OO 

k 

855 

*  0 .  31433E+00 

0 . 28570E+00 

0  -OOOOOE+OO 

k 

856 

*  0.39891E+00 

0 . 14841E+00 

0. OOOOOE+OO 

k 

857 

*  -0 . 18211E+00 

0 . 38615E+00 

0. OOOOOE+OO 

k 

858 

*  0.20978E+00 

0  -37325B+00 

0.  OOOOOE+OO 

k 

859 

*  -0 .  40851E+00 

0.13171E+00 

0. OOOOOE+OO 
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*  S60  *  0 . 39812E+00  0.1651QE+00  O.OOOOOE+OG 

*  361  *  0 . 28632E+00  3.32250E+00  0.00000E+00 

*  362  *  -  0 . 32385E+00  0.28570E+00  O.OOOOOE+OO 

*  863  *  0 . 39540E  +  00  0.18225E+00  0.00000E  +  00 

’  364  *  -  0 . 29584E+00  0.32250E+00  O.OOOOOE+OO 

'  365  *  -  0 . 22883E  +  00  0.37325E  +  00  O.OOOOOE  +  OO 

*  366  *  0 . 25471E+00  0.35707E+00  O.OOOOOE+OO 

*  367  *  -  0 . 40702E  +  00  0.16510E.00  O.OOOOOE  +  OO 

'  368  *  0 . 39201E  +  00  0.19926E+00  O.OOOOOE  +  OO 

’  369  *  O.OOOOOE+OO  0.44331E+00  O.OOOOOE+OO 

*  870  *  0 . 12452E-01  0.44327E+00  O.OOOOOE+OO 

*  871  *  -  0 . 12452E -  01  0.44327E  +  00  O.OOOOOE+OO 

*  872  *  -  0 . 40429E+00  0.18225E+00  O.OOOOOE+OO 

*  873  *  0 . 24904E -  01  0.44315E+00  O.OOOOOE+OO 

*  874  *  -0 . 24904E-01  0.44315E+00  O.OOOOOE+OO 

*  875  *  - 0 . 26424E+00  0.35707E+00  O.OOOOOE+OO 

*  876  *  0 . 37356E- 01  0.44295E+00  O.OOOOOE+OO 

*  877  *  -  0 . 37356E- 01  0.44295E+00  O.OOOOOE+OO 

*  878  *  0 . 44493E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  379  *  0 . 49808E-01  0.44266E+00  O.OOOOOE+OO 

*  880  *  -0 . 49808E-01  0.44266E+00  O.OOOOOE+OO 

*  881  *  0 . 38415E+00  0.22760E+00  O.OOOOOE+OO 

*  882  *  0 . 62260E- 01  0.44230E+00  O.OOOOOE+OO 

*  883  *  -  0 . 62260E- 01  0.44230E+00  O.OOOOOE+OO 

*  884  *  0.44451E+00  0.49150E-01  O.OOOOOE+OO 

*  885  *  -  0 . 40090E+00  0.19926E+00  O.OOOOOE+OO 

*  886  *  0 . 74711E- 01  0.44186E+00  O.OOOOOE+OO 

*  887  *  - 0 . 74711E- 01  0.44186E+00  O.OOOOOE+OO 

*  888  *  0 . 87163E-01  0.44134E+00  O.OOOOOE+OO 

*  889  *  -0 . 87163E-01  0.44134E+00  O.OOOOOE+OO 

*  890  *  0 . 99615E- 01  0.44074E+00  O.OOOOOE+OO 

*  891  *  -0.99615E-01  0.44074E+00  O.OOOOOE+OO 

*  892  *  0 . 37413E+00  0.25515E+00  O.OOOOOE+OO 

*  893  *  -  0 . 45319E+00  O.OOOOOE+OO  O.OOOOOE+OO 

*  894  *  -0 . 11437E+00  0.43897E+00  O.OOOOOE+OO 

*  895  *  0 . 44324E+00  0.98290E-01  O.OOOOOE+OO 

*  896  *  0 . 12081E+00  0.43780E+00  O.OOOOOE+OO 

*  897  *  - 0 . 39304E+00  0.22760E+00  O.OOOOOE+OO 

*  898  *  - 0 . 12904E+00  0.43675E+00  O.OOOOOE+OO 

*  899  *  -0 .45276E+00  0.49148E-01  O.OOOOOE+OO 

*  900  *  0 . 14179E+00  0.43375E+00  O.OOOOOE+OO 

*  901  *  -0.14115E+00  0 . 43503E+00  O.OOOOOE+OO 

*  902  *  0 . 36204E+00  0.28059E+00  O.OOOOOE+OO 

*  903  *  0 . 16383E+00  0.42862E+00  O.OOOOOE+OO 

*  904  *  -  0 . 15321E+00  0.43309E+00  O.OOOOOE+OO 

*  905  *  - 0 . 38302E+00  0.25515E+00  O.OOOOOE+OO 

*  906  *  0 . 44188E+00  0.13172E+00  O.OOOOOE+OO 

*  907  *  0 . 18551E+00  0.42243E+00  O.OOOOOE+OO 

*  908  *  -0.45149E+00  0.98290E-01  O.OOOOOE+OO 

*  909  *  -0 . 17396E+00  0.42822E+00  O.OOOOOE+OO 

*  910  *  0 . 34798E+00  0.30489E+00  O.OOOOOE+OO 

*  911  *  0 . 20680E+00  0.41520E+00  O.OOOOOE+OO 

*  912  *  -  0 . 19440E+00  0.42243E+00  O.OOOOOE+OO 

*  913  *  - 0 . 37093E+00  0.28059E+00  O.OOOOOE+OO 

*  914  *  0 . 33204E+00  0.32725E+00  O.OOOOOE+OO 

*  915  *  0 . 22760E+00  0.40695E+00  O.OOOOOE+OO 

*  916  *  0 . 31437E+00  0.34814E+00  O.OOOOOE+OO 

*  917  *  - 0 . 35687E+00  0.30489E+00  O.OOOOOE+OO 

*  918  *  -0.22024E+00  0.41520E+00  O.OOOOOE+OO 

*  919  *  0.44014E+00  0.16510E+00  O.OOOOOE+OO 

*  920  *  0 . 25243E-!  00  0.39772E+00  O.OOOOOE+OO 

*  921  *  - 0 . 34093E+00  0.32725E+00  O.OOOOOE+OO 
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0 . 29607E+00 

0 .36863E+00 

923 

Hr 
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0 . 34814E+00 
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Hr 
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0 .40695E+00 
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Hr 
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0 . 38754E+00 
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Hr 

-0 .44839E+00 
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★ 

-  0 . 26566E+00 

0 .39772E+00 
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Hr 

-  0 . 30496E+00 

0 .36863E+00 
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★ 

0 . 43301E+00 

0 .2Q352E+00 

930 

Hr 

-0 .28534E+00 

0 .38754E+00 

931 

Hr 

0. OOOOOE+OO 

0 . 48307E+00 

932 

Hr 

0 . 25456E- 01 

0 . 48288E+00 

933 

Hr 

-  0 . 25456E -  01 

0 . 48288E+00 

934 

Hr 

0 . 50912E-01 

0 . 48232E+00 

935 

Hr 

-  0 . 50912E- 01 

0 . 48232E+00 

936 

Hr 

-  0 . 44126E+00 

0 . 20352E+00 

937 

H* 

0 . 76368E- 01 

0 . 48139E+00 

938 

Hr 

-0 . 76368E-01 

0 . 48139E+00 

939 

Hr 

0 .48763E+00 

0 .OOOOOE+OO 

940 

Hr 

0 .48751E+00 

0 . 24576E-01 

941 

★ 

0 . 48715E+00 

0 . 49150E-01 

942 

Hr 

0 . 10182E+00 

0 . 48007E+00 

943 

★ 

-  0 . 10182E+00 

0 . 48007E+00 

944 

Hr 

0 . 41214E+00 

0 . 26711E+00 

945 

Hr 

0 .48654E+00 

0 . 73721E-01 

946 

Hr 

-  0 . 13468E+00 

0.47542E+00 

947 

Hr 

0 . 14955E+00 

0 .47192E+00 

948 

★ 

-0.49525E+00 

0. OOOOOE+OO 

949 

★ 

0 . 48569E+00 

0.98290E-01 

950 

★ 

-0 . 49513E+00 

0 . 24575E- 01 

951 

★ 

-0 .49477E+00 

0.49148E-01 

952 

Hr 

-0 . 16075E+00 

0 . 47115E+00 

953 

Hr 

-0.42039E+00 

0 . 26711E+00 

954 

Hr 

0 . 48497E+00 

0 . 11501E+00 

955 

★ 

-0.49416E+00 

0.73720E-01 

956 

Hr 

0 . 19844E+00 

0.45872E+00 

957 

Hr 

0 .38164E+00 

0 . 32408E+00 

958 

H- 

0  .48414E+00 

0 . 13172E+00 

959 

Hr 

-0.49331E+00 

0 . 98290E-01 

960 

* 

-0.20670E+00 

0.45872E+00 

961 

★ 

0 . 24543E+00 

0 -44066E+00 

962 

★ 

0 . 48320E+00 

0 . 14841E+00 

963 

* 

0 . 34242E+00 

0 . 37378E+00 

964 

Hr 

-  0 . 38989E+00 

0.32408E+00 

965 

★ 

-0.49176E+00 

0 . 13172E+0Q 

966 

★ 

0.48215E+00 

0 . 1651QE+00 

967 

Hr 

-0 . 35067E+00 

0 . 37378E+00 

968 

★ 

-  0 . 26194E+00 

0.44066E+00 

969 

* 

0.29818E+00 

0.41801E+00 

970 

★ 

0 . 47852E+00 

0.18653E+00 

971 

Hr 

-0.48977E+00 

0.16510E+00 

972 

Hr 

0 . 47400E+00 

0.20777E+00 

973 

Hr 

-0.30643E+00 

0 . 41801E+00 

974 

Hr 

-  0 . 48614E+00 

0 . 18653E+00 

975 

★ 

0.  OOOOOB+OO 

0 . 52283E+00 

976 

* 

0.13004B-01 

0 . 52278E+00 

977 

* 

-0 . 13004E-01 

0.52278E+00 
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* 

0.26008E-01 

0.52262E+00 

979 

Hr 

-0.26008E-01 

0 . 52262E+00 
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Hr 
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981 

★ 
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0 . 52235E+00 

982 

Hr 

-0.39012B-01 
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983 

Hr 
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984 

* 

0 .52016E-01 

0  -52198E+00 

985 

* 

-  0 . 52016E- 01 

0 .52198E+00 

986 

* 

0 . 65021E -  01 

0 . 52150E+00 

987 

★ 

-0 . 65021E-01 

0 .52150E+00 

988 

★ 

0 . 78025E -  01 

0 . 52091E+00 

989 

* 

-0 . 78025E-01 

0 .52091E+00 

990 

* 

0 . 91029E- 01 

0 . 52021E+00 

991 

+ 

-0 .91029E-01 

0 . 52021E+00 

992 

•k 

0 .45015E+00 

0 . 27907E+00 

993 

* 

0 . 10403E+00 

0 .51941E+00 

994 

★ 

-  0 . 10403E+00 

0 . 51941E+00 

995 

★ 

0 . 53033E+00 

0. OOOOOE+OO 

996 

*■ 

-0 .47113E+00 

0 . 24395E+00 

997 

★ 

-  0 . 12223E+00 

0 .51706E+00 

998 

-* 

0 . 13081E+00 

0 .51549E+00 

999 

* 

0 .52978E+00 

0 .49150E-01 

1000 

★ 

-0 .14031E+00 

0 .51409E+00 

1001 

■k 

0 . 15730E+00 

0 . 51009E+00 

1002 

* 

0 . 43403E+00 

0 . 31194E+00 

1003 

■k 

-  0 . 15434E+00 

0 . 51179E+00 

1004 

0 . 18457E+00 

0 . 50326E+00 

1005 

★ 

-  0 . 45777E+00 

0 .27907E+00 

1006 

* 

-0 .16830E+00 

0 .50921E+00 

1007 

*• 

0 .52815E+00 

0 .98290E-01 

1008 

★ 

-0 .53731E+00 

0 .OOOOOE+OO 

1009 

* 

0 . 21137E+00 

0 .49500E+00 

1010 

•k 

-  0 . 19384E+00 

0 . 50271E+00 

1011 

* 

0 . 41529E+00 

0 .34328E+00 

1012 

■k 

-0 .53677E+00 

0 . 49148E-01 

1013 

* 

0 . 23763E+00 

0.48536E+00 

1014 

★ 

-0 .44165E+00 

0 . 31194E+0Q 

1015 

★ 

-0 .21899E+00 

0 .49500E+QQ 

1016 

★ 

0 .39403E+00 

0 . 37233E+00 

1017 

★ 

0 . 26325E+00 

0 . 47436E+00 

1018 

* 

0 . 52641E+00 

0 . 13172E+00 

1019 

* 

-  0 . 53513E+00 

0 . 98290E-01 

1020 

* 

-  0 . 42291E+00 

0 . 34327E+00 

1021 

it 

0 . 37047E+00 

0 . 39942E+00 

1022 

* 

-0 .24915E+00 

0 . 48536E+00 

1023 

★ 

0 . 29207E+00 

0 . 46205E+00 

1024 

- 0 . 40165E+00 

0 . 37233E+00 

1025 

★ 

0 . 34608E+00 

0 . 42500E+00 

1026 

★ 

0 . 52416E+00 

0 . 16510E+00 

1027 

* 

-  0 . 37809E+00 

0.39942E+00 

1028 

★ 

-  0 . 27849E+00 

0 .47436E+00 

1029 

★ 

0 . 31991E+00 

0.44848E+00 

1030 

* 

-0 . 30341E+00 

0.46205E+00 

1031 

★ 

-0.35370E+00 

0 . 42500E+00 

1032 

* 

-  0 . 32753E+00 

0.44848E+00 

1033 

★ 

-  0 . 53115E+00 

0 . 16510E+00 

1034 

★ 

0 . 51500E+00 

0 . 21203E+00 

1035 

★ 

0. OOOOOE+OO 

0 . 56260E+00 

1036 

# 

0 . 26560E-01 

0 . 56236E+00 

1037 

★ 

-0.26560E-01 

0 . 56236E+00 

1038 

★ 

-0 . 52198E+00 

0 .21203E+00 

1039 

★ 

0 . 53121E-01 

0 . 56164E+00 

1040 

* 

-0.53121E-01 

0 . 56164E+00 

1041 

★ 

0 . 79681E-01 

0 . 56043E+00 

1042 

★ 

-  0 . 79681E-01 

0 . 56043E+00 

1043 

★ 

0 . 48817E+00 

0 . 29103E+00 

1044 

★ 

0 . 10624E+00 

0 . 55875E+00 

1045 

-0 .10624E+00 

0 . 55875E+00 

0 . OOOOOE+OO 
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G .OOOOOE+OO 
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0.57302E+00 
0 . 57287E+00 

- 0 . 49515E+00 
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*  1108  *  - 0 . 18338E+00 

*  1109  *  0 . 23722E+G0 

'  1110  *  0 . 4826QE-rOG 

*  1111  *  -0.21373E+00 

*  1112  *  0 . 61S72E+00 

*  1113  *  -  0  .  S1237E  +  00 

*  1114  *  0 . 26846E+00 

*  1115  '  0 . 61506E+00 

*  1116  *  -  0 . 2435  8E  +  00 

*  1117  *  0 . 45603E+00 

*  1118  *  0 . 29890E+00 

*  1119  *  -  0 . 48895E+00 

*  1120  *  0 . 42658E+00 

’  1121  *  0 . 61306E+00 

*  1122  *  - 0 . 27806E+G0 

*  1123  *  -  0 . 62144E+00 

*  1124  *  0 . 33171E+00 

*  1125  *  - 0 . 62077E+00 

*  1126  *  -  0 . 46238E+00 

*  1127  *  0 . 39609E+00 

*  1128  *  - 0 . 43292E+00 

*  1129  *  0 . 61093E+00 

*  1130  *  -  0 . 31161E+00 

*  1131  *  0 . 36338E+00 

*  1132  *  -  0 . 61878E+00 

*  1133  *  - 0 . 34116E+00 

*  1134  *  -  0 . 40244E+00 

*  1135  *  -0.3  6972E+00 
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*  1140  *  0  .  Q00Q0E+0Q 
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*  1147  *  0 . 56419E+00 

*  1148  *  0 . 11066E+00 

*  1149  *  -  0 , 11Q66E+00 

*  1150  *  -0.1S721E+00 
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+  1642  *  - 0 . 99728E+00  0.73717E-01  O.OOOOOE+OO 

*  1643  *  0 . 99879E+00  0.49150E-01  O.OOOOOE+OO 

*  1644  *  -  0 . 99879E+00  0.49145E-01  O.OOOOOE+OO 

*  1645  *  0 . 99970E+00  0.24576E-01  O.OOOOOE+OO 

*  1646  *  -0.99970E+00  0.24573E-01  O.OOOOOE+OO 

*  1647  *  - 0 . 10000E+01  O.OOOOOE+OO  O.OOOOOE+OO 

*  1648  *  0 . 1000QE+01  O.OOOOOE+OO  O.OOOOOE+OO 

*  1649  *  -0.98893E+00  0.14843E+00  O.OOOOOE+OO 

*  1650  *  - 0 . 99336E+00  0.11503E+00  O.OOOOOE+OO 
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677 

643 

588 

685 

633 

*369* 

693 

793 

677 

769 

744 

685 

783 

728 

*370* 

793 

892 

769 

868 

845 

783 

881 

823 

*371* 

892 

992 

868 

972 

944 

881 

980 

929 

*372* 

992 

1096 

972 

1077 

1043 

980 

1079 

1034 

*373* 

1096 

1200 

1077 

1181 

1147 

1079 

1183 

1138 

*374* 

1200 

1302 

1181 

1285 

1251 

1183 

1287 

1242 

*375* 

1302 

1406 

1285 

1389 

1355 

1287 

1391 

1346 

*376* 

1406 

1510 

1389 

1493 

1459 

1391 

1495 

1450 

*377* 

1510 

1626 

1493 

1630 

1561 

1495 

1628 

1554 

*378* 

583 

677 

573 

655 

633 

578 

665 

622 

*379* 

677 

769 

655 

756 

728 

665 

759 

711 

*380* 

769 

868 

756 

860 

823 

759 

863 

812 

*381* 

868 

972 

860 

966 

929 

963 

970 

919 

91 


*382* 

972 

1077 

*383* 

1077 

1181 

*384* 

1181 

1285 

*395* 

1285 

1389 

*386* 

1389 

1493 

*387* 

1493 

1630 

*388* 

573 

655 

*389* 

655 

”56 

*  39C* 

756 

860 

*391* 

860 

966 

*392* 

966 

1072 

*393* 

1072 

1177 

*394* 

1177 

1281 

*395* 

1281 

1385 

*396* 

1385 

1489 

*397* 

1489 

1634 

*398* 

541 

641 

*399* 

641 

747 

*400* 

747 

850 

*401* 

850 

958 

*402* 

958 

1066 

*403* 

1 0  6  o 

1173 

*404* 

1173 

1278 

*405* 

1278 

1381 

*406* 

1381 

1485 

*407* 

1485 

1637 

*408* 

510 

624 

*409* 

624 

736 

*410* 

736 

841 

*411* 

841 

949 

*412* 

949 

1056 

*413* 

1056 

1166 

*414* 

1166 

1272 

*415* 

1272 

1376 

*416* 

1376 

1480 

*417* 

1480 

1640 

*418* 

464 

598 

*419* 

598 

715 

*420* 

715 

831 

*421* 

831 

941 

*422* 

941 

1051 

*423* 

1051 

1162 

*424* 

1162 

1270 

*425* 

1270 

1374 

*426* 

1374 

1478 

*427* 

1478 

1643 

*428* 

404 

414 

*429* 

414 

469 

*430* 

383 

400 

*431* 

400 

447 

*432* 

360 

379 

*433* 

373 

433 

*434* 

352 

369 

*435* 

369 

418 

*436* 

351 

368 

*437* 

368 

417 

*438* 

335 

348 

*439* 

348 

398 

*440* 

322 

337 

*441* 

337 

381 

*442* 

315 

329 

*443* 

329 

2  71 

966 

1072 

1034 

1072 

1177 

1138 

1177 

12  81 

1242 

1281 

1385 

1346 

1385 

1489 

1450 

1489 

1634 

1554 

541 

641 

622 

641 

”47 

”11 

747 

850 

312 

350 

958 

919 

958 

1066 

1026 

1066 

1173 

1136 

1173 

1278 

1240 

1278 

1381 

1344 

1381 

1485 

1448 

1485 

1637 

1552 

510 

624 

587 

624 

736 

692 

736 

841 

799 

841 

949 

906 

949 

1056 

1018 

1056 

1166 

1129 

1166 

1272 

1238 

1272 

1376 

1342 

1376 

1480 

1446 

1490 

1640 

1550 

4  64 

598 

567 

598 

715 

674 

715 

831 

784 

831 

941 

895 

941 

1051 

1007 

1051 

1162 

1121 

1162 

1270 

1235 

1270 

1374 

1339 

1374 

1478 

1443 

1478 

1643 

1547 

448 

591 

539 

591 

708 

653 

708 

824 

761 

824 

939 

884 

939 

1048 

999 

1048 

1159 

1115 

1159 

1268 

1233 

1268 

1372 

1338 

1372 

1476 

1442 

1476 

1648 

1546 

448 

464 

408 

464 

510 

439 

404 

414 

385 

4X4 

469 

416 

383 

400 

365 

400 

447 

402 

360 

379 

354 

379 

433 

391 

335 

348 

353 

348 

398 

390 

322 

337 

339 

337 

381 

362 

315 

329 

325 

329 

371 

350 

311 

323 

320 

323 

363 

346 

970 

1075 

1026 

1075 

1179 

1136 

1179 

12  8  3 

1240 

1283 

1387 

1344 

1387 

1491 

1448 

1491 

1632 

1552 

559 

646 

587 

646 

”50 

692 

750 

356 

799 

356 

962 

906 

962 

1068 

1018 

1068 

1176 

1129 

1176 

1280 

1238 

1280 

1383 

1342 

1383 

1487 

1446 

1487 

1635 

1550 

523 

632 

567 

632 

739 

674 

739 

346 

784 

346 

954 

895 

954 

1062 

1007 

1062 

1170 

1121 

1170 

1275 

1235 

1275 

1372 

1339 

1379 

1482 

1443 

1482 

1638 

1547 

483 

606 

539 

606 

721 

653 

721 

835 

761 

835 

945 

884 

945 

1053 

999 

1053 

1163 

1115 

1163 

1271 

1233 

1271 

1375 

1338 

1375 

1479 

1442 

1479 

1641 

1546 

452 

593 

532 

593 

709 

647 

709 

825 

760 

825 

940 

878 

940 

1049 

995 

1049 

1161 

1112 

1161 

1269 

1230 

1269 

1373 

1337 

1373 

1477 

1441 

1477 

1645 

1545 

427 

441 

452 

441 

489 

483 

393 

412 

408 

412 

460 

439 

373 

389 

385 

389 

437 

416 

358 

375 

365 

375 

429 

402 

344 

356 

339 

356 

410 

362 

331 

341 

325 

341 

387 

350 

318 

333 

320 

333 

377 

346 

313 

327 

316 

327 

367 

342 
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’444  * 

311 

323 

314 

328 

316 

312 

226 

319 

’445* 

323 

363 

328 

270 

342 

326 

366 

345 

*446* 

314 

328 

321 

336 

319 

317 

332 

324 

’447* 

328 

370 

236 

380 

345 

3  32 

376 

349 

*448* 

321 

336 

234 

24  7 

324 

330 

340 

338 

-449* 

336 

380 

347 

397 

349 

340 

386 

361 

*450* 

334 

247 

352 

3  69 

338 

343 

355 

354 

*451* 

347 

397 

269 

418 

361 

355 

409 

391 

*452* 

359 

378 

351 

368 

364 

357 

374 

353 

*453* 

378 

432 

368 

417 

401 

374 

428 

390 

*454  * 

382 

399 

359 

378 

3  84 

372 

388 

364 

’455* 

399 

446 

378 

432 

415 

388 

436 

401 

*456* 

403 

413 

382 

399 

407 

292 

4  11 

334 

*457* 

413 

468 

399 

446 

438 

411 

459 

415 

*458* 

449 

465 

403 

412 

453 

426 

440 

407 

*459* 

465 

511 

413 

468 

484 

440 

438 

438 

*460* 

474 

490 

449 

465 

478 

463 

477 

453 

*461* 

490 

527 

465 

511 

509 

477 

518 

484 

*462* 

1613 

1544 

1617 

1539 

■>  570 

1£1  r 

1543 

1576 

463* 

1544 

1440 

1539 

1435 

1475 

1543 

1439 

1472 

*464* 

1440 

1336 

1435 

1331 

1371 

1439 

1335 

1368 

*465* 

1336 

1234 

1331 

1227 

1267 

1335 

1232 

1264 

*466* 

1234 

1135 

1227 

1128 

1174 

1232 

1134 

1164 

*467* 

1135 

1032 

1128 

1027 

1074 

1134 

1031 

1069 

*468* 

1032 

930 

1027 

923 

973 

1031 

928 

967 

*469* 

930 

834 

923 

819 

875 

928 

828 

864 

*470* 

834 

737 

819 

722 

782 

828 

731 

762 

*471* 

737 

634 

722 

623 

681 

731 

630 

668 

*472* 

1617 

1539 

1621 

1533 

1576 

1619 

1537 

1574 

*473* 

1539 

1435 

1533 

1429 

1472 

1537 

1433 

1470 

*474* 

1435 

1331 

1429 

1325 

1368 

1433 

1329 

1366 

*475* 

1331 

1227 

1325 

1221 

1264 

1329 

1225 

1262 

*476* 

1227 

1128 

1221 

1119 

1164 

1225 

1126 

1158 

*477* 

1128 

1027 

1119 

1020 

1069 

1126 

1024 

1063 

*478* 

1027 

923 

1020 

917 

967 

1024 

921 

964 

*479* 

923 

819 

917 

814 

864 

921 

818 

862 

*480* 

819 

722 

814 

714 

762 

818 

717 

758 

*481* 

722 

623 

714 

617 

668 

717 

621 

656 

*482* 

1621 

1533 

1625 

1520 

1574 

1623 

1528 

1567 

*483* 

1533 

1429 

1520 

1417 

1470 

1528 

1424 

1463 

*484* 

1429 

1325 

1417 

1313 

1366 

1424 

1320 

1360 

*485* 

1325 

1221 

1313 

1210 

1262 

1320 

1216 

1256 

*486* 

1221 

1119 

1210 

1106 

1158 

1216 

1113 

1152 

*487* 

1119 

1020 

1106 

1005 

1063 

1113 

1014 

1050 

*488* 

1020 

917 

1005 

905 

964 

1014 

913 

953 

*489* 

917 

814 

905 

804 

862 

913 

810 

853 

*490* 

814 

714 

804 

702 

758 

810 

707 

751 

*491* 

714 

617 

702 

605 

656 

707 

613 

652 

*492* 

1625 

1520 

1629 

1496 

1567 

1627 

1511 

1555 

*493* 

1520 

1417 

1496 

1392 

1463 

1511 

1407 

1451 

*494* 

1417 

1313 

1392 

1288 

1360 

1407 

1305 

1347 

*495* 

1313 

1210 

1288 

1184 

1256 

1305 

1201 

1243 

*496* 

1210 

1106 

1184 

1083 

1152 

1201 

1097 

1139 

*497* 

1106 

1005 

1083 

983 

1050 

1097 

996 

1038 

*498* 

1005 

905 

983 

885 

953 

996 

897 

936 

*499* 

905 

804 

885 

788 

853 

897 

796 

838 

*500* 

804 

702 

788 

689 

751 

796 

697 

742 

*501* 

702 

605 

689 

597 

652 

697 

602 

642 

*502* 

1629 

1496 

1633 

1492 

1555 

1631 

1494 

1553 

*503* 

1496 

1392 

1492 

1388 

1451 

1494 

1390 

1449 

*504* 

1392 

1288 

1388 

1284 

1347 

1390 

1286 

1345 

*505* 

1288 

1184 

1284 

1180 

1243 

1286 

1182 

1241 

93 


*506* 

1184 

1083 

1180 

1076 

1139 

1182 

1078 

1137 

*507* 

1083 

983 

1076 

971 

1038 

1078 

974 

1033 

*508* 

983 

385 

971 

367 

936 

574 

372 

926 

*509* 

385 

■788 

867 

”66 

838 

372 

-78 

321 

*510* 

7S8 

689 

•766 

5  78 

”4  2 

””8 

682 

”29 

*511* 

689 

597 

578 

586 

542 

682 

592 

635 

*512* 

1633 

1492 

1639 

1488 

1553 

16  36 

1490 

1551 

-513* 

1492 

1388 

1488 

1384 

1449 

1490 

1386 

1447 

-514* 

1388 

1284 

1384 

1279 

1345 

1386 

1282 

1343 

*515* 

1284 

1180 

1279 

1175 

1241 

1282 

1178 

1239 

*516* 

1180 

1076 

1175 

1067 

1137 

1178 

1073 

1132 

*517* 

1076 

971 

1067 

959 

1033 

1073 

965 

1019 

*518* 

971 

867 

959 

852 

926 

965 

859 

908 

*519* 

967 

766 

852 

745 

821 

859 

753 

798 

*520* 

766 

678 

745 

638 

729 

753 

649 

688 

*521* 

678 

586 

638 

527 

635 

649 

555 

577 

*522* 

1639 

1488 

1644 

1484 

1551 

1642 

1486 

1549 

*523* 

1488 

1384 

1484 

1380 

1447 

1486 

1382 

1445 

*524* 

1384 

1279 

1380 

1276 

1343 

1382 

1277 

1341 

*525* 

1279 

1175 

1276 

1171 

1239 

1277 

1172 

1237 

*526* 

1175 

1067 

1171 

1061 

1132 

1172 

1065 

1125 

*527* 

1067 

959 

1061 

951 

1019 

1065 

955 

1012 

*528* 

959 

852 

951 

844 

908 

555 

348 

399 

*529* 

852 

745 

844 

733 

798 

348 

738 

787 

*530* 

745 

638 

733 

618 

688 

738 

629 

671 

*531* 

638 

527 

618 

490 

577 

629 

509 

556 

*532* 

1644 

1484 

1647 

1481 

1549 

1646 

1483 

1548 

*533* 

1484 

1380 

1481 

1377 

1445 

1483 

1378 

1444 

*534* 

1380 

1276 

1377 

1273 

1341 

1378 

1274 

1340 

*535* 

1276 

1171 

1273 

1168 

1237 

1274 

1169 

1236 

*536* 

1171 

1061 

1168 

1057 

1125 

1169 

1059 

1123 

*537* 

1061 

951 

1057 

948 

1012 

1059 

950 

1008 

*538* 

951 

844 

948 

839 

899 

950 

840 

893 

*539* 

844 

733 

839 

725 

787 

840 

730 

779 

*540* 

733 

618 

725 

607 

671 

730 

609 

660 

*541* 

618 

490 

607 

474 

556 

609 

478 

547 

LINE03R  FLUIDE  3 

*542*  1648  1643  1645 
*543*  1643  1640  1641 
*544*  1640  1637  1638 
*545*  1637  1634  1635 
*546*  1634  1630  1632 
*547*  1630  1626  1628 
*548*  1626  1622  1624 
*549*  1622  1618  1620 
*550*  1618  1614  1616 
*551*  1614  1609  1611 
*552*  1609  1605  1607 
*553*  1605  1601  1603 
*554*  1601  1595  1598 
*555*  1595  1591  1593 
*556*  1591  1587  1589 
*557*  1587  1583  1585 
*558*  1583  1580  1581 
*559*  1580  1584  1582 
*560*  1584  1588  1586 
*561*  1588  1592  1590 
*562*  1592  1596  1594 
*563*  1596  1599  1597 
*564*  1599  1602  1600 
*565*  1602  1606  1604 
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*566* 

1606 

1610 

1608 

*567* 

1610 

1613 

1612 

*568* 

1613 

1617 

1615 

*569* 

1617 

1621 

1619 

'570* 

1621 

1625 

1623 

*571* 

1625 

1629 

1627 

*572* 

1629 

1633 

1631 

*573* 

163  j 

1636 

1619 

*574* 

1636 

1639 

1650 

*575  * 

1639 

1644 

1642 

*576* 

1644 

1647 

1646 

END 

“  i 

3 

3 

-1 

4 

2 

s 

'  3 

-4 

-2 

*  Non  existing  degree  of  freedom  ;.2-D) 

*  Electrode  at  0  V 

*  Electrode  at  i  V 
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APPENDIX  B 


MATLAB  PROGRAM  FOR  ANALYTICAL  MODEL  OF  UNLOADED  TRANSDUCER 


%  This  MATLAB  program  was  writcen  by  MAJ  Tay  Tiong  Beng 
%  on  X  Feb  1993  to  implement  the  formula  developed 
%  by  Kong- zhang  Wang  [Ref.  2]  for  unloaded  tangentially 
s'  polarized  ring  transducer 
% 

%  Initialization 

clear ; 
clg ; 

!del  res 
!del  resl.met 
:del  res2.met 
!del  res3.met 
diary  res 

%  Material  Constants 

j  =sqrt ( - 1)  ; 
rhoc=7002 ; 

sll=0 . 1702e- 10*  ( 1 -0 . 015* j )  ;  %  Losses  is  1.5%  of  the  real  part.  The  real  part 

%  is  the  published  value  of  the  elastic  constant 

sl2=-0.405e-ll; 

nu=- (sl2) / (sll)  %  Poisson  Ratio 

s66=2* (S11-S12) 

e33=0 . 57e- 8 

d31=-Q ,135e-9 

d33=0 . 3e-9 

Yc=l/sll ; 

F=3000 ;  %  Frequency 

%  Fluid  Constants 

cair=343 ; 
cwater=1500 ; 
ccer=sqrt (Yc/rhoc) ; 

%  Dimensions 

a=  (14 . 44/2) *0 . 0254  %  Inside  radius 

b=  (14 .94/2) *0.0254  %  outside  radius  of  ceramic  stave 

1=6.8*0.0254  %  height 

%  Capacitance  and  Admittance 

N=72  ; 

Co=N*2*e33*l/ (2*pi) *log(b/a) 
k31s=Yc*(d31*2) /e33 
t=(d33/d31) 

Cs=Co* (1- (k31s/ (l-nu*2) ) * (l+2*nu*t+tA2) ) 


V  Admittance  Circle 


3 =sqrt  t - 1) ; 

wr=ccer/'a  %  Radial  Res 

fr=wr / ( 2  *pi i 
n  =  0  ; 

for  J  =  1 : 5  :  F 
n=n+l ; 
w=2*pi*f ; 
wwsw/wr ; 

WW ( 1 ,  n) =ww ; 
j wcs=j *w*Cs ; 

fa= tl-nuA2) * !w/wr) *2  -  1; 

A= j *w«Co*k3 Is / f a ; 

ka2= (w/ccer) "2*fa/ ( (w/wr) A2-l5 ; 

KA2  (1 ,  n>  =ka2  ; 
ka=sqrt (ka2 ) ; 

KA ( 1 , n) =ka; 

B= (fa* (l+nu*t) +nu* (nu+t) ) A2/ ( (l-nu'2) A2* ( (w/wr) A2-l) )  ; 

D= (nu+t) "2/ (l-nuA2) ; 

T=A*  (B* (tan (0 . 5*ka*l) / (0 ,5*ka*l! ) -D) +  :jwcs; 

TT=T/ ( j  *wr*Co)  ; 

2=1/ TT ; 

Z  ( 1 ,  n )  =  z  ; 

Y(l,n) =T ; 

YY  ( 1 ,  n)  =TT  ; 

end; 

q=abs ( Y) 

G=real (Y) 

S=imag (Y) 

f =1 : 5 : F ; 

axis ( ' square' )  ; 

plot (G, S) ; grid; 

title (’ADMITTANCE  CIRCLE  OF  A  TANGENTIALLY  POLARIZED  UNLOADED  TRANSDUCER') 
xlabel (' Conductance ,  G  mhos'); 
y label ( ' Susceptance ,  S  mhos' ) ; 
meta  resl ; 

plot (f ' ,G) ;  grid; 

title ( 'CONDUCTANCE  AS  A  FUNCTION  OF  FREQUENCY'); 

xlabel (' Frequency.  Hz');  ylabel (' Conduct ance ,  G  mhos');  pause; 

meta  res2 ,- 

plot ( f ' , S) ;  grid; 

title (' SUSCEPTANCE  AS  A  FUNCTION  OF  FREQUENCY' ) ; 

xlabel (' Frequency,  Hz');  ylabel (' Susceptance,  S  mhos');  pause; 

meta  res3 ; 

diary  off 
end; 
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APPENDIX  C 


TYPICAL  VALUES  FOR  MATERIALS  CONSTANTS 


TABLE  Is  PHYSICAL  DIMENSIONS  AND  MATERIAL  CONSTANTS 


The  Transducer 

Inside  Radius  :  0.1651  m 

Outside  Radius  :  0.204724  m 

Height  :  0.2921  m 

Motor  Element:  Navy  Type  1  Ceramic  (Hexagonal  Class  6mm) 

Density  :  7550  kg/m3 

Inside  Radius  :  0.183388  m 

Outside  Radius  :  0.189738  m 

Elastic  Constants  (m2/N) 

S33e  =  1.539E-11 
s,3e  -  -5.31E-12 
Si2e  =  -4.05E-12 
S„E  =  1.22E-11 
S44E  »  3.90E-11 
s66E  =  3.25E-11 

Dielectric  Constants 

633t  =  1.149E-8 
€  i  jT  =  1.2963E-8 
633s  =  0.57E-8 
eus  =  0.624E-8 


Piezoelectric  Constants 
g  tensor  (m2/c) 

g33  =  2  6.  IE -3 
g3l  =  -11.7E-3 
g15  =  40.5E-3 

d  tensor  (m/v) 

d„  =  300E-12 
d31  »  -135E-12 
dj5  -  525E-12 
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TABLE  1:  CONTINUED 


Fiberglass  Wrap 


Thickness 

Density 

Longitudinal  Modulus 
Transverse  Modulus 
Poisson  Ratio 


0.762  mm 
2  540  kg/m3 

9 . 5E6  psi  ; 5 . 8 6 E 1 0  Pa) 
2.9E6  psi  (2.00E10  Fa) 
0.28 


Encapsulant  Material 


Polyurethane 

Density  :  0.039  pci  {  >=1100  kg/nr) 

Young's  Modulus  :  600  psi  (4.14  MPa) 

Poisson  Ratio  :  0.49 

Neoprene 

Density  :  0.045  pci  (  »1250  kg/m3) 

Young's  Modulus  :  100-3000  psi 

(0.69-20.69  MPA) 

Poisson  Ratio  :  0.49 

Reflector  Rina:  HRS  Carbon  Steel 

Density  :  7700  Kg/m3 

Young's  Modulus  :  19.5E10  Pa 
Poisson  Ratio  :  0.28 

Mounting  Plate:  Fiberglass -epoxy  composite 

Density  :  1700  Kg/m3 

Young  Modulus  :  2000  psi  (13.79  MPa) 

Poisson  Ratio  :  0.34 
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APPENDIX  D 


SAMPLE  CALCULATIONS  FOR  THE  EFFECTIVE 
ELASTIC  CONSTANTS  OF  THE  FIBERGLASS  WRAPPED  CERAMIC 

RING 


The  effective  elastic  constants  of  the  ring  ire  estimated  as 
follows.  Consider  a  section  of  the  fiberglass  wrapped  ring  as  two 
rectangular  volume  elements  cemented  together  as  shown  m  Figure  1 
below.  The  purpose  is  to  calculate  the  elastic  constants  of  an 
equivalent  ceramic  material  which  will  give  the  same  elastic 
properties  as  the  fiberglass  wrapped  ceramic  ring.  The  equivalent 
ceramic  material  has  the  same  length,  1,  and  height,  h,  as  the 
fiberglass  and  ceramic,  and  a  width  equal  to  the  sum  of  the  widths 
of  the  fiberglass  and  ceramic. 


Figure  1 
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The  thickness  of  the  fiberglass  is  t  and  that  of  the  ceramic  i 
where  t  and  w  are  0. 000762m  and  0.0064m,  respectively. 

Consider  first  an  applied  force,  F},  in  the  circumferen 
direction  only,  as  shown.  Assuming  that  the  resultant  stress 
is  uniform  then  the  strain,  5,,  in  each  material  is  given  by: 


For  the  ceramic, 


(Ah),.  =  h.(sJ3)t.  Ts 


And  for  the  fiberglass, 


(Ah)?  =  hf(s13)  f  T} 


There  is  no  reason  to  assume  (Ah)  g  equals  to  (Ah)c.  We 
therefore  define  (Ah) eff  for  the  equivalent  material  as  the  average 
Ah  across  the  composite  material,  neglecting  boundary  effects: 

(Ah)  cw+ (Ah)  g  t  =  (A h)eff  (w+t)  (4) 


Then  the  effective  strain  in  the  axial  direction,  i.e.,  in  the  "1" 
direction  is  given  by: 


(S)  -  {Ah)*  - 

x)«  ~h~ 


(S13)cw+(sl3)tt 


w+  t 


after  taking  hc  ~  ha  =  h. 
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Thus  the  effective  compliance, 


snj  erf •  is  giver,  by: 


i 


13'  r’X 


sl}l  w-  s:J. ,  t 
_____ 


(6) 


By  considering  forces  m  the  ocher  directions,  in  turn,  the 
same  result  is  obtained,  with  the  appropriate  identification  of 
indices  and  element  dimensions.  For  example,  to  compute  s,:  we  will 
consider  an  applied  force  in  the  ’'2”  direction,  i.e.,  the  radial 
direction,  and  calculate  the  strain  in  the  " 1 " ,  or  axial, 
direction. 


Using  values  listed  in  Appendix  C  for  the  fiberglass  and 
ceramic,  the  following  were  obtained: 


CONSTANTS 

. 

COMPUTED  VALUES  (m*/N) 

St  |  =  S->2 

1.142  x  10'" 

S\2  *  S21 

-4.9  x  10‘i2 

S,3  =  S31 

-5.3  x  10i2 

S23  ~  Si2 

-5.3  x  10‘12 

S33 

1.702  x  10  “ 

I5 

II 

U* 

3.9  x  10‘u 

^66 

3.58  x  10*n 
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